i 


B 


TE 


ee 


Na 
me 


=. 


Digitized by the Internet Archive 
in 2023 with funding from 
Kahle/Austin Foundation 


https://archive.org/details/paulilecturesonp0001 paul 


: Pauli Lectures on Physics: 
_ Electrodynamics 


s*+2 


Pauli Lectures on Physics 


1. Electrodynamics 

2. Optics and the Theory of Electrons 

3. Thermodynamics and the Kinetic Theory of Gases 
4, Statistical Mechanics 

5. Wave Mechanics 


6. Selected Topics in Field Quantization 


“Cheese 


1 2 
at 


I | Pauli Lectures on Physics: 
' Volume 1. Electrodynamics 


Wolfgang Pauli 


Edited by Charles P. Enz 
Translated by S. Margulies and H. R. Lewis 


Foreword by Victor F. Weisskopf 


Library 
Western Wyoming Community Golegs 


The MIT Press 
Cambridge, Massachusetts, and London, England 


Copyright © 1973 by 
The Massachusetts Institute of Technology 


All rights reserved. No part of this book may be reproduced in any form 
or by any means, electronic or mechanical, including photocopying, 
recording, or by any information storage and retrieval system, without 
permission in writing from the publisher. 


First MIT Press paperback edition, 1977 
Second printing, 1978 
Third printing, 1979 


ISBN 0 262 66033 4 (paperback) 
Library of Congress catalog card number: 76—155320 


iis o 


Contents 


Foreword by Victor F. Weisskopf 
Preface by the Editor 


Survey of the Historical Development and the Current 


Problems of Electrodynamics 
1. Electrostatics and Magnetostatics 


1. Coulomb’s law 

2. The field of point charges 

3. Volume and surface charge distributions 

4. Energy of the electrostatic field 

5. Example: charge distributions with spherical symmetry 
6. Proof of the equivalence of the electrostatic field equations 

with Coulomb’s law 

7. Dielectrics, phenomenological treatment 

8. Electron theory interpretation of dielectric phenomena 
9. The potential problem 
10. Curvilinear coordinates 

11. Examples of solutions of the potential problem 

12. Magnetostatics 

13. Units and dimensions 


Steady-State Currents 


14. Theory of steady-state currents 
15. The magnetic fields of steady-state currents 
16. The equivalence of the magnetic field due to closed current 
loops with that due to a distribution of magnetic dipoles 
Vv 


vii 


1x 


vi ; CONTENTS 


17. Ponderomotive forces 77 
18. The principle of action and reaction for electric and 

magnetic forces. The Maxwell stress tensors 81 
3. Quasi-Static Fields 88 - 
19. Faraday’s law of induction 88 
20. The energy of current systems 90 
21. Time-dependent current flow in circuits 94 
22. The skin effect 103 
23. The law of induction for moving conductors 106 
4. Rapidly Varying Fields 109 
24. The Maxwell equations 109 
25. Electromagnetic waves in vacuum 113 
26. Conservation of energy and momentum 117 
27, Electromagnetic waves in material media 124 
28. Radiation of electromagnetic waves 125 
29. The field of a point charge in uniform motion 143 
30. Radiation damping 146 
Bibliography 153 
Appendix. Comments by the Editor 155 


om 


Index 157 


Foreword 


It is often said that scientific texts quickly become obsolete. 
Why are the Pauli lectures brought to the public today, when 
some of them were given as long as twenty years ago? The 
reason is simple: Pauli’s way of presenting physics is never 
out of date. His famous article on the foundations of quantum 
mechanics appeared in 1933 in the German encyclopedia 
Handbuch der Physik. Twenty-five years later it reappeared 
practically unchanged in a new edition, whereas most other 
contributions to this encyclopedia had to be completely re- 
written. The reason for this remarkable fact lies in Pauli’s 
style, which is commensurate to the greatness of its subject 
in its clarity and impact. Style in scientific writing is a quality 
that today is on the point of vanishing. The pressure of fast 
publication is so great that people rush into print with hur- 
riedly written papers and books that show little concern for 
careful formulation of ideas. Mathematical and instrumental 
techniques have become complicated and difficult; today most 
of the effort of writing and learning is devoted to the acquisi- 
tion of these techniques instead of insight into important 
concepts. Essential ideas of physics are often lost in the dense 
forest of mathematical reasoning. This situation need not be 
so. Pauli’s lectures show how physical ideas can be presented 
clearly and in good mathematical form, without being hidden 
in formalistic expertise. 

Pauli was not an accomplished lecturer in the technical sense 
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of the word. It was often difficult to follow his courses. But 
when the sequence of his thoughts and the structure of his 
logic become apparent, the attentive follower is left with a new 
and deeper knowledge of essential concepts and with a clearer 
insight into the splendid architecture of reason, which is theo- 
retical physics. The value of the lecture notes is not diminished 
by the fact that they were written not by him but by some of 
his collaborators. They bear the mark of the master in their 
conceptual structure and their mathematical rigidity. Only 
here and there does one miss words and comments of the mas- 
ter. Neither does one notice the passing of time in his lectures, 
with the sole exception of the lectures on field quantization, 
in which some concepts are formulated in a way that may 
appear old-fashioned to some today. But even these lectures 
should be of use to modern students because of their compact- 
ness and their direct approach to the central problems. 

May this volume serve as an example of how the concepts of 
theoretical physics were conceived and taught by one of the 
great men who created them. 


Victor F. Weisskopf 


Cambridge, Massachusetts 


Preface 


As stated in the introductory survey, Pauli chose to teach 
electrodynamics in the inductive way, in opposition to his 
teacher, Arnold Sommerfeld, who had preference for the axi- 
omatic presentation of Maxwell’s equations (see the preface 
in A. Sommerfeld, Electrodynamics, Academic Press, New 
York, 1952). This attitude is typical of Pauli and can also be 
found in the other lectures of this series (with the exception 
of field quantization). It expresses Pauli’s vivid interest in the 
formation of scientific concepts and ideas and the logical struc- 
tures built upon them—an eminently historic process (see, 
e.g., the introduction in W. Pauli, “Der Einfluss archetypischer 
Vorstellungen auf die Bildung naturwissenschaftlicher The- 
orien bei Kepler” in Naturerklarung und Psyche, Rascher 
Verlag, Ziirich, 1952). 

Because of this concern for the logical structure of the 
theory, the present lectures are still a rewarding source of 
learning for students of today. The notes prepared by A. Thel- 
lung, on which this English translation is based, reflect well 
the concise style of Pauli’s lecturing. When these notes were 
published at ETH, Ziirich, in 1949, Thellung was Pauli’s Ph.D. 
student. He later became Pauli’s assistant and is now Professor 
of theoretical physics at Zurich University. The care and pre- 
cision of his notes, so typical of Thellung, made the work of 
the translators a comparatively easy job. 

For Pauli the central problem of electrodynamics was the 
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field concept and the existence of an elementary charge which 
is expressible by the fine-structure constant e?/me = 1/137. 
This fundamental pure number had greatly fascinated Pauli, 
as can be seen from the list of references to his work assembled 
in the appendix. For Pauli the explanation of the number 137 
was the test of a successful field theory, a test which no theory 
has passed up to now. This number 137 transcended into a 
magic symbol at Pauli’s death. When I visited Pauli in the 
hospital, he asked me with concern whether I had noticed his 
room number: 137! It is in this room that he died a few days 
later. 


Charles P. Enz 


Geneva, 17 November 1971 


Survey of the Historical Development and the 
Current Problems of Electrodynamics 


Electrodynamics is a relatively young branch of theoretical 
physics. As a domain of field physies it has become very 
important in recent times. Field physics is concerned with 
continuous functions in space and time (examples: elasticity 
theory, hydrodynamics, etc.), in contrast to corpuscular 
physics, which also plays a role in electrodynamics since 
electricity has an atomistic structure. 

The concept of a field goes back to Faraday, who employed 
an intuitive approach without rigorous mathematical for- 
mulas, and who introduced the concept of lines of force. 
Maxwell brought the theory into a rigorous mathematical 
form. The existence of the electromagnetic waves predicted 
by this theory was later verified experimentally by Hertz. 
Since light is nothing but electromagnetic waves within a 
certain wavelength region, it became possible to consider 
optics as a branch of electrodynamics. 

The mechanics of continuous media provided the model 
for field physics. Thus, by analogy, it was believed that an 
“aether” was required for the propagation of light, and an 
attempt was made to attribute the phenomena of the prop- 
agation of light to the mechanical properties of this aether. 
(Mechanistic concept of nature: all the laws of nature were 
to be explained on the basis of mechanics.) Maxwell, Boltz- 
mann, and others worked out so-called aether engines. How- 
ever, it was soon seen that all these theories became very arti- 
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ficial and complicated. Thus, little by little, this approach 
was abandoned and the concept of a physical field became 
firmly established. However, electromagnetic fields can be 
measured only with so-called test charges; hence they are 
linked with concepts of corpuscular physics. 

It has become evident that the aether cannot possess any 
atomistic structure (as indeed there is in the case of elastic 
media), and that the concept of motion can be applied only 
to the test charge and not to the aether itself. This latter 
fact has led to the development of the theory of relativity. 

It is, however, by no means true that field physics has 
triumphed over corpuscular physics, as can be seen from 
the fact that electricity is atomistic in nature. The carriers 
of electrical charge are negative electrons and positive pro- 
tons. In addition, positive electrons (positrons), negative 
antiprotons, and other particles (mesons and hyperons) of 
both positive and negative charge can be artificially pro- 
duced by suitable processes. All electrical charges are in- 
tegra] multiples of an elementary charge. 

H. A. Lorentz and J. Larmor traced the differences in 
the electrical behavior of macroscopic bodies, characterized 
by so-called ‘‘material constants,’ back to Maxwell’s equa- 
tions in vacuum and to carriers of electric charge (electron 
theory). They thus achieved a great simplification of the 
foundations, which proved to be of special advantage for rel- 
ativity theory. On the other hand, there is no explanation 
for the fact that only integral multiples of a certain charge 
occur. The existence of an elementary charge has, until 
now, in no way been made plausible. It is still an open 
problem in theoretical physics. The electron itself is a stran- 
ger in the Maxwell-Lorentz theory as well as in present-day 
quantum theory [A-1]. 2 

The field-particle description presents a conceptual prob- 
lem: although a field can be described mathematically 


‘Comments [A-1]-[A-4] appear in the Appendix on pp. 155 -156. 
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without the need for any test charges, it cannot be measured 
without them. On the other hand, the test charge itself 
gives rise to a field. However, it is impossible to measure 
an external field with a test charge and, at the same time, 
to determine the field due to this charge. A certain duality 
exists. Consequently, electrodynamics is of great significance 
for physical epistemology [A-1]. 

Electrodynamics can be presented in two ways: 
1. Deductive: starting with the Maxwell equations and 
developing special cases. 
2. Inductive: beginning with the basic laws obtained from 
experiment and eventually building up to the Maxwell 
equations. This way corresponds more closely to the 
historical development. 
In these lectures, we will employ the second approach. 


Chapter 1. Electrostatics and Magnetostatics 


1. COULOMB’S LAW 


The concept of conductors and insulators is assumed to 
be familiar from experimental physics. In these paragraphs 
we will consider only charged bodies whose linear dimen- 
sions are small as compared with their relative separation. 
Thus, we can assume that their electric charges are con- 
centrated at points (point charges). 

The forces between two such point charges are central 
and satisfy the principle of action and reaction (action=re- 
action). The forces can be either attractive or repulsive 
and their intensities are given by 


C= |e, °€, | 


ye? 


where e, and e, are the electric charges and r is their sep- 
aration. In vector notation (Fig. 1.1), the force on the 


€2 
This 
ey 
x: 
X 
O 


Figure 1.1 
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charge e, due to the charge e, is 
Cama 
ae aaa [1.1] 


where x, and x, are the position vectors of e, and e, 

Ty2=|*,—x,| is the separation between e, and @, my 

=(%_—4%,)/r,. is a unit vector in the direction from e, to e. 

This empirically obtained relation is called Coulomb’s law. 
One must add two other experimentally deduced facts: 

1. Superposition of forces. Electric forces add according to 
the laws of vector addition (parallelogram rule): 


Ky +03 = K, 3+ K, ;. 


2. Conservation of charge. Electric charges, which can be 
neither created nor destroyed, add algebraically. 

Then Coulomb’s law implicitly provides a sufficient de- 

finition of electrical charge since various quantities can 

be varied. 

For example, one defines: (a) two charges to be equal 
when each produces the same effect upon a third charge; 
(b) a charge e, to be twice as large as charge e, if it pro- 
duces the same effect upon a third charge e, as do two 
charges, each of magnitude e,. The distance between these 
two charges e, must, of course, be small as compared with 
their distance to e,. Electric charges can be either posi- 
tive or negative. Charges of like sign repel, those of op- 
posite sign attract. 

One can also consider the forces exerted by a definite 
charge e, on various charges e,, and vice versa. This dual 
aspect forms the basis for the definition of the electric field: 
For a given e¢,, the force on a test charge e, is 


Kish [1.2] 


Here E, the electric field intensity, is independent of ¢. 
For the case of a field originating from a single point 
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charge eg at the point Q, the field at a point P is given 
by Coulomb’s law: 


é 6, oo6(1 
E, =— == (f= 2 — 6. Ls 
iP t Nop aes aie @) Q A) Xp oe b [ ] 


where Q, the source point, is fixed, and P, the field point, 
is arbitrary. For the case of several point charges, it fol- 
lows from the principle of superposition that 


€ € 
E,= ¥-2 nop = i ye 1.4 

. 2, a 22, canis Oxp ie L1-4) 
In the general case, for any arbitrary charge distribution, 
the electric field intensity is also defined by the equation 


K=cE. [1.5] 


Here K is the force on the test body, e¢ is the electric charge 
carried by the test body and depends only on the test 
body, and E is the electric field intensity which is inde- 
pendent of the test body and depends only on the field. 

It is essential that we assume the existence of test charges 
that do not perturb the field or whose effects on the field 
are negligibly small. In macroscopic electrostatics this pro- 
duces no difficulties since e, the magnitude of the test 
charge, can be made practically arbitrarily small. Further- 
more, in the case that the sources of the field are fixed 
point or quasi-point charges, e can even be arbitrarily large. 
On the other hand, if the sources of the field are extended 
charged conducting surfaces, then a large test charge would 
produce a displacement of the charges, thereby changing 
the field (see See. 11). Contrary to macroscopic electro- 
statics, great difficulties arise in the atomic domain. Here, 
the reaction of the test charge on the field cannot be ne- 
glected since the test charge cannot be made arbitrarily 
small and, moreover, since the field sources are not at rest. 
This points to a definite difficulty in the field concept [A-1]. 
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Up to this point we have described the field in the lan- 
guage of corpuscular physics. In the following sections we 
will consider the field defined by Coulomb’s law as an 
independent concept and will investigate its properties. 


2. THE FIELD OF POINT CHARGES 


a. The field is conservative (conservation of energy) 
From vector analysis we know that this property can be 
expressed by any of four fully equivalent expressions: 


i Eds 0 [2.1] 
for an arbitrary closed curve. The physical significance is 


eh Eds =p K-ds—0, 


that is, no net work is either gained or lost when a charged 
particle traverses a closed path in an electric field. 


ie 
2. [E-as is dependent only upon the position of the end 

o 
points O and P (Fig. 2.1) and not on the path between 


tg 


Figure 2.1 


them. From this it can easily be shown that this integral 
must have the form 


[Bas=—g9 +90. [2.2] 
o 
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Proof: We define (E -ds= F(O, P). Then, for an arbitrary 
point P’, ‘ 
F(0, P’) + F(P', P) = F(O, P) 
a FO, PP) — FP, Po =TOeey, 
since obviously F(P’,P)=—F(P, P’). We now keep P’ 
fixed: 
F(O, P')=q, and F(P, P')=¢,. 
Therefore, 
Yo— Pr = FO, P), 


which is the desired result. 

The quantity gy, is called the electrostatic potential at the 
point P. It is defined only to within an arbitrary additive 
constant. (The zero of potential can be chosen arbitrarily.) 


3. E = — grado. [2.3] 
In other notations, 
Ee VG; ie a=W Zt, t= 1.2 
4, curl — 0) . [2.4] 
Other notations are 
Of; 0o&#; 


VxE=0; 


aa Ga i, j cyclic. 
] cy 


The equivalence of condition 4 with the three other forms 
of the energy law follows from Stokes’s theorem: 


[Eas = curl E-n df. [2.5] 


For this equation, in a right-handed coordinate system, the 
relation between the normal direction and the sense of 
rotation of C shown in Fig. 2.2 must be maintained. 
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From each of the above four statements the other three 


follow. They can be shown to be valid for any arbitrary 
central force. Thus, they are also valid for a finite sum 


n 
df 
C F 
ds 
Figure 2.2 


of central forces and, consequently, also for the field of 
point charges. 
The potential due to a system of point charges is 


[2.6] 


Here, we have chosen the arbitrary additive constant so 
that g(co) =0. The correctness of this expression for g can 
be shown by forming the gradient: 
E = — gradp9 = > Sf (xp— x0) 
@ Top 
(The differentiation is with respect to the coordinates of P; 
the sources remain fixed.) 

Since statements 1, 2, 3, and 4 are valid for arbitrary 
potential fields, they say far less about the electric field E 
than does Coulomb’s law. They are equivalent to Coulomb’s 
law only when taken together with the following law. 


b. Electric field flux. Gauss’s law 


1. A single point charge. The electric field flux through a 
sphere with the source point at its center (Fig. 2.3) is 


f #. af =| 5rd = Ane. 
r2 
x x 
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The flux is the same for all spheres regardless of their size. 
One can easily see by extension that for an arbitrary closed 
surface 


f Bua =—h aj 


PF 


__|4ne if @ lies inside the surface F, - 
~ |0 if Q lies outside the surface F. [2.7] 
K 

nm 
Figure 2.3 


2. Several point charges. In'this case, 
) 
p Bndf=— PE af = tn Sey. [2.8] 
n @ 
F F 


The sum is to be-taken over all the charges contained within 
the surface F. 
Because of Gauss’s theorem, 


pi, af= p En aj =[aivEar, [2.9] 


in which F is a closed surface which bounds the volume V 
and n is the outward normal, it follows from Eqs. [2.7] 
and [2.8] that 


divpE=0 for P+Qq. [2.10] 
Since E = —gradq@ and 


Cee et ae 
div grad » = V’o = Ay= (sat ant eal ' 
3 
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g satisfies Laplace’s equation: 
Vep=0 for P#Q. [2.11] 


It is immediately seen that y=e/r (for a single point 
charge) or y= > ¢,/r,, (for several point charges) are solu- 
@ 


tions of this equation, because 


Cr hy 0 ()= L; oO (7) ==" 


5 = ==> —-|]— 
On, On; \7r re? Oar 


(1) & fl\_ 3r—3r 
vf) Pel) a= 


except at the origin. For a surface containing the origin, 


Pali sari m5) ay = 4a. 


Energy conservation, 


curlE=0 or E=-~—gradg, 


and Gauss’s law, 


p Bu df = 42 Dey or ~ 2 aj = 4230, 
@ on @ 


F 


(Q within F) are, together, completely equivalent to 
Coulomb’s law for point charges. 


Mathematical remarks on Stokes’s and Gauss’s theorems 


In the case of Stokes’s theorem we had a definite relation 
between the direction of the normal n and the sense of tra- 
versal of the curve C. This convention is necessary only if 
curl E and the surface element df are considered as vectors. 
It is, however, more consistent to consider them as anti- 
symmetric tensors.! This can be seen by considering an 


1 For the concept of a tensor, see Section 18. 
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inversion of the coordinate system. However, to any anti- 
symmetric tensor in a three-dimensional space, 


Cu = A,B,— A,B; SS Cis 


(the A; and B, are to transform as vector components), a 
vector C, can be assigned with the help of the tensor é:,,, 
which is antisymmetric in all three indices: 


+1 if tkl is an ea permutation of 1,2,3, 


Ext = 


0 if at least two of the three indices are equal. 
That is, one sets 


C;= > EniCrr. 
k<l 

Then, the vector components C,, C,, and C,;are equal to the 
tensor components C,;, C3, and C,,, respectively, in all 
coordinate systems which can be derived by proper rotations 
(i.e., where the determinant of the rotation matrix is +1) 
of the original system. If improper rotations (with deter- 
minant —1) are considered, then the signs change. For 
example, 


~ 


f 


Ll; =— 2; (¢=1, 2,3), 
On= Crs but 0; _ Oe . 


From the point of view of inversions, both the curl and 
the surface element can be considered as antisymmetric 
tensors (vector products). Thus, one can write 


pE-as = » curl,, Edf i. . [2.12] 
t<k 

Hence, an ordering between the normal direction and the 

sense of rotation is not necessary. Indeed, the concept of 

a normal never appears. In this form Stokes’s theorem is 
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valid in an arbitrary n-dimensional space. This very general 
theorem is purely topological; that is, it is independent of 
any metric. It is only necessary that the H, transform con- 
tragrediently to the da’. That is, 


E-ds = > E, da‘ = invariant. 


A closed surface can always be thought of as being com- 
posed of two separate surfaces with a common boundary. 
If Eq. [2.12] is applied to both of these surfaces, then the 
line integrals along the bounding curves cancel and it fol- 
lows, for a closed surface, that 


> Meurl, E df, = 0. 
t<k 
The proof of Stokes’s theorem proceeds most easily by 
first proving it for a rectangle and then showing that the 
theorem is invariant under arbitrary coordinate transfor- 
mations. Thus, it must be true for all surfaces that can be 
transformed into a rectangle. For the case of multiply- 
connected regions, a separation into singly-connected regions, 
for each of which the theorem then applies, is made by 
appropriate cuts. Since the line integrals cancel in pairs 
along the cuts, the theorem is thus also valid for multiply- 
connected regions. 
Starting with a vector A; and an antisymmetric tensor 
Cx=—Cyi, one can construct the tensor 


A; Cat A; Cin+ A, Cu = Dini . 


Since D,,, is antisymmetric in all three indices it can always 
be written in the form 


Dia = Fim D « 


In Gauss’s theorem, the concept of the normal component 
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can now be eliminated by writing 
p (E,Ofi+ Ey fis + Bi dfx) 
F 


= | div E(AV ) sx: =|divE | AV | ésx1 . [2.13] 
V 


It is seen that Gauss’s theorem is also of topological char- 
acter. The left-hand side transforms as a determinant, as 
does the right-hand side, since divE is a scalar. As is true 
for Stokes’s theorem, Gauss’s theorem is tied neither to a 
3-dimensional space nor to the concept of a metric. 

As in the case of Stokes’s theorem, Gauss’s theorem is 
most easily proved by showing it to be true for a cube, 
then showing that it is invariant under arbitrary coordi- 
nate transformations, and so forth. 


3. VOLUME AND SURFACE CHARGE DISTRIBUTIONS 


It is often more convenient to think of charge distributions 
as being continuous instead of consisting of a series of point 
charges. As long as one remains in the domain of macro- 
scopic physics, the actual atomic structure of electricity 
can be neglected. 


a. Volume charges 


For the charges contained within a volume V we imagine 
the following limiting process: 


the number of charges — oo, 


the magnitude of each individual charge +0, 


to take place in such a manner that the total charge remains 
constant. Thus, 


2 €—g—>e={0dV (total charge in Y). 


V 
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The continuous function thus defined, o= lim (e/ V), is called 


the volume charge density. The potential g, of Eq. [2.6] and 
the field E, of Eq. [1.4] then become 


af fear, in 


E, = — grad> 9p =| £ (Xp— XQ) dg. [3.2] 


In both of these integrals there are no singularities as the 
point P approaches Q, contrary to what occurs in the sum- 
mation of point charges. Stated very roughly, the integrands 
OQo/"op and 0,(x,—x,)/7>, become infinite more slowly than 
the volume element dV approaches zero. This can easily be 
shown to be rigorously true by an exact consideration of 
the limit. Thus, for continuous volume distributions of 
charge, E and 9 are regular inside the source region. 
For the case of volume charges, Gauss’s law becomes 


) 
—f 2 af=o Haj = safeav [3.3] 


Because of the regularity of the integrand, this can, with 
the aid of Gauss’s theorem (Eq. [2.9]), be transformed into 


-[v'pav—|aivEav— tn| oa. 
V v v 


Since this is true for any arbitrary volume V, and in par- 
ticular for arbitrary small volumes, then 


div E = 470 [3.4] 
or 

— Vp = 4n0. [3.5] 
Equation [3.5], which is a generalization of Laplace’s equa- 
tion (Eq. [2.11]), is called Poisson’s equation. 
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b. Surface charges 


One can imagine a process exactly analogous to that 
employed for volume charges applied to a two-dimensional 
surface. Thus, one arrives at the concept of a surface charge 
density w. In analogy, 


e=|wdf (total charge on the surface) , 


PF 
(43) 
— | aah, [3.6] 
QP 
(39) 
E, =— gradpy =| (xp— XQ) dia - [3.7] 


As the field point P passes into the surface, y remains finite 
(integrand ~1/r, surface element ~r?). Thus p is continuous 
at a charged surface. 

On the other hand, this is not the case for the compo- 
nent of the electric field normal to the surface. This is seen 
by applying Gauss’s law of Eq. [2.9] to an imaginary 


df: _n, 


Nee df, 


Figure 3,1 


cylindrical surface containing a section of the charged 
surface, a8 shown in Fig. 3.1: 


pk, af= taf w af . [3.8] 


In the limit, as the height of the cylinder goes to zero while 
the base remains constant, 


feet Bay = taf wa 
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Since the base area is arbitrary, 
E,, + E,, = 40. : [3.9] 


Alternately, with a different definition of the normal 
(Fig. 3.2), 
(En): — (En)z = 4 , 


Op Op 
(=2|— (=), = — 47w. 


Thus, the normal component of E is discontinuous at a 
charged surface. On the other hand, the tangential com- 


Figure 3.2 


ponent remains continuous. Proof: y is the same on either 
side of the surface. Thus, differentiation parallel to the 
surface in the same direction must yield the same result, 
both inside and out. The proof can also be made with 
Stokes’s theorem. 


Boundary conditions for the field at conductors: Electrical 
conductors are substances with freely movable electrons. 
Electric foree components inside a conductor produce 
charge displacements—a current flows. The conditions for 
no current flow or, equivalently, for the static character 
of the field, are 
1. no field may be present inside the conductor: 


= 0; 


2. at the surface of the conductor the tangential component 
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of the field vanishes: 
Ei — 0 . 


(Since the tangential component is continuous, it must be 
zero outside the conductor just as it is inside.) On the con- 
trary, the normal component outside the conductor is, in 
general, not zero since surface charges may be present. 
Thus, E inside is zero while externally the field consists of 
just the normal component 


Be Ee ay. 
on 


The gradient of the potential is perpendicular to the 
surface. Thus, g =constant along this surface. That is, 
the surface is an equipotential surface. 


4. ENERGY OF THE ELECTROSTATIC FIELD 


a. Point charges 


The potential energy of a system of point charges is equal 
to the work required to bring the charges from infinity to 
their actual configuration. For two charges e, and e, at a 
distance rj, the work is 


BE, =22 


2 
pot = 
12 


For an arbitrary number of point charges, by calculating the 
potential energy for each pair and summing over all pairs, 


€,e 
Ew = > —. [4.1] 
i<k Vk 
One can also write 
1 €,é 
BE eel ¢-%ivtk ' 
pot = 5 > re [4.2] 


The prime on the summation symbol is to indicate that 
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terms with i= are not included. From this symmetrical 
expression it can be seen that the energy thus defined is 
independent of the order in which the point charges are 
brought from infinity. 

There is still a third form of £,,, in which the concept of 
potential is used. The potential at the position of the ith 


charge due to all charges except the ith is 


y= =. [4.3] 
k ik 

(k#%) 
Consequently, 

all 


b. Volume charges 


ee 5 / | Grr ay ay, . [4.5] 
Tpp’ 


Here, the complication of i#k, which exists for point 
charges, does not appear, since in the expression [ (Open) V, 
it makes no difference whether the volume element contain- 
ing P is included or not. 

In the notation using the potential, 


p= | ave, [4.6] 


Ba = 5 [enedVe. [4.7] 


In this case 


c. Volume charges and surface charges 


With the potential 


Qe -| Oa: dV +f aed df» (Eqs. [4.6] plus [3.6]), 
Tp’ Ypp' 


there results 
1 1 
Eesot = 5 er Qp adVp+ 3 WpYp df, . [4.8] 
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The first integral extends over all regions which contain 
volume charge distributions, the second over all charged 
surfaces. Without an explicit introduction of the potential 
concept, 


Eyot a L Cree dV, dv’, +f = dV, df, 
2 Ypp’ Y pp’ 


alt Wp Wp 
He | Or dfn dfe. [4.9] 


As a consequence of the laws of electrostatics, these expres- 
sions for energy can be transformed into integrals over all 
space. 

First, for volume charges only, using Eq. [3.4], 


Moot 5 erav— rag divE dV. 


1 
Now, with the familiar identity from vector analysis, 
div (pE) = pdivE + E-gradg, [4.10] 


which is valid for an arbitrary vector E and an arbitrary 
scalar function y, Gauss’s theorem yields 


pos a pdivEay +[E-gradgay : 


provided that div E exists over the whole region of integ- 
ration. (Thus, E must be continuous.) Applied to our sit- 
uation, this results in 


il 
Ego = a {-[#-eraapar + p ob. aj] : 
,.¢ 


As the region of integration we imagine a sphere K whose 
radius 7 approaches infinity in the limit. We now assume 
that @~#0 only within a finite sphere, or else that @ van- 
ishes sufficiently fast at infinity so that the second integral 
over the sphere K vanishes as K->co. (If @40 only 
within a finite region, then m0 at least as fast as 1/r 
and |E|—0 at least as 1/r?, while the surface increases 
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only as 7°. Thus, the whole integral approaches zero at 
least as 1/r.) Because E=—gradg, one thus obtains 


il 
Biya = = | EVV , [4.11] 


in which the energy is expressed in terms of only the field in- 
tensity at all points in space. The quantity E?/8x is called 
the energy density of the field. 

If surface charges are also present, the corresponding sur- 
faces F must be excluded from the region of integration since 


Figure 4.1 


there E is not continuous and Gauss’s theorem is not appli- 
cable. This gives us two additional terms which, however, 
cancel: 


1 il 
Ba = 3 [eva +5 coy af 
1 
= {- | B-eraag) av + pebadl + pb aa +5 fowar. 
7-4 F 
Now, 


il 1 =f 

es a + E,,) df =— = |oedf, 

sap PBoa a |. Jaf Oy alll j 
F 


since the normal n, which must point away from the region 
of integration, is oppositely directed to m, and n, (see 
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Fig. 4.1). If we again assume that the integral over the 
sphere vanishes as K->oo, then 


1) 8 fe dv [4.11] 
87 


is again obtained. This formula is thus valid for both volume 
and surface charges. 

On the other hand, this integra] diverges for point charges. 
In this case, the field must be decomposed into components 
arising from the various point charges. The potential and 
the field at the point P produced by the point charge e, are 

pr P) = = » E,=—gradg,. 
We now consider the integral [E:-E,.aV for i~k. If we 
exclude the singularities by enclosing the point charges 


Uj 
GG 


GG 


Figure 4.2 


within small spheres K; and K,, then, remembering that 
div E=0 over the region of integration, Gauss’s theorem 
yields 


[EE av= ~[E.-gradp, dV 


- ~[aiv (Bay +| .divE,dV 


7 ~ f Bapeat — h Bosgea f Ba Px Of . 
E Ex 


Ki 
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The first term approaches zero as the large sphere becomes 
infinite. The third term approaches zero as we let K, go 
to zero. Only the second term contributes. For K;—0, 


_— p Ei, af — Ane; . 
xy 


Thus, 


a $ Ein Qr df = An Caer 


Vix 
For two point charges, therefore, 


67,0; 
Take 4a 


oot => E;:E, av * [4.12] 


Here, the integral extends over all space. 
For an arbitrary number of point charges 


ey eo? __ E;- -E,AV= = > | E,-E,dV. [4.13] 
=72 ik 


i<k Vix 


One can also write 


eee | SENav. [4.14] 


(In this case, however, the integral cannot be decomposed 
into two individual integrals since these would diverge.) 


5. EXAMPLE: CHARGE DISTRIBUTIONS WITH SPHERICAL SYMMETRY 


We consider a volume distribution of charge o(r) which 
depends only upon the distance r from a fixed point 0. 
Thus, the potential g(r) is a function of r only. 
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In this case it is simplest to start from Poisson’s equation 
— V°p = 4x0. [5.1] 


In polar coordinates the operator V?, when applied to a 
function g dependent only upon r, has the form 


V2 =s5(" 2-7" as ae” [5.2] 


dr} rdr 


d 
—nt = tal or dr =e(r), 
0 
where e(7) is the charge contained within a sphere of radius r. 
Hence, 
ap _ par 


a = [5.3] 


The field at a point located a distance r from O, the center 
of the charge distribution, depends only upon the charge 
contained within a sphere of radius r about O. The charge 
outside does not affect E as long as it is distributed with 
spherical symmetry. The field is the same as it would be 
if all of the charge e(r) were concentrated at the center O. 
The potential is given by 


vr) =—[ Dar = ena (=) 
i —[Fewar. 


= e(r) i 


Because e'(r) = (d/dr)[e(r)] =4ar?e, and since, for r—0, 
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e(r) >0 as r?, we have 
= 


G(r) = o — te grdr+ constant . [5.4] 


0 


a. Uniformly charged spherical volume 


Oo = constant for r<a 
v= 5 
0 fOr 7 >a; 
4 rs 
Bets — oo, r<a 
e(r) = 
(r)=) 4 
te Oye Cea, 
va 
ap a (ats } 
Sadr OTT eg 
P r>a, 
é 
qa (8a — 7") Waa: 
i 
e 
- roa 


(The constants have been chosen so that y(co) =0 and » 
is continuous at r=a.) Calculation of the electrostatic 


energy yields 


1 NH NG 
Evo = 5 fer a 3 5 5 | (eat ryder dr 
0 


: aft 
The same result is obtained by using = Er av. 
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b. Uniformly charged spherical surface 
Let 
_ @ 
0 ara? 


where ¢ is the total charge and a istheradius of the sphere. 
As a consequence of Gauss’s law, 


ry? B= constant 


must be valid both inside and outside the sphere. If we 
denote the region inside the sphere by 1 and the region 
outside by 2, then 


Cy 
r<a: E,=—, 
7? 


but c,= 0, since there can be no discontinuity at the center. 
Hence 


Ey, = 0 

€ 
r>a: EB, = —; 
v2 


but 


é é 

r>a Et Ges) 
é 
r<a: Ee= 0, o=-=. 
a 


Again, we see that the surface charge acts externally as 
though it were concentrated at the center of the sphere. 
Within the sphere it produces no effect. The energy is 


For the limiting case of a point charge, E,,,—> oo. 
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6. PROOF OF THE EQUIVALENCE OF THE ELECTROSTATIC FIELD 
EQUATIONS WITH COULOMB’S LAW 


We have derived the field equations [2.4], [2.3], or [3.4] 

[3.5], 
curl E = 0, div E = 4x, 
or 
E=— grad, V2p = — 400, 

from Coulomb’s law. We now wish to show that, con- 
versely, Coulomb’s law follows from these equations. In 
order to demonstrate this we must, in addition, require 
that y—0 at least as fast as 1/r as roo. 

For simplicity we will perform the proof only for vol- 
ume charges. For this, we will need Green’s theorem. If 
A= ggrady is substituted into Gauss’s theorem, 


pAndf =p A-naf = div AdV, 
F F Lie 


then, because 
div (p grad y) = gp V*p + grady-grady, [6.1] 


we obtain Green’s first identity: 


post al— pV'pdV + [erady-grady aV. [6.2] 


Interchanging @ and y in [6.1] we obtain, upon subtraction, 
Green’s second identity: 


div (py grad yp — p grad) = p Vp — p V29, 
) rs) 
b(o5¢- 4 | df =[(eV~— Vp) AV, [6.3] 


These results are valid for any arbitrary volume in which 
gy and y are regular. 

Green’s second identity can be applied to the integration 
of Poisson’s equation by considering g as the sought po- 


28 : ELECTROSTATICS AND MAGNETOSTATICS | Chap. 1 


tential function and letting yp =1/rpp. 1n order that the 
integrands be regular, we limit our integration to the region 


P' (variable) 


Figure 6.1 


between a small sphere Kp about the fixed point P and a 
large sphere K (see Fig. 6.1). Because 


V4yv=—0> for eer’, 
we have 


O(i/r) 1 
fyi 12 


a ae aa 
+ (9 ane OR et 


KP 


Now, the integral over the sphere K vanishes because of 

our assumption about the behavior of » at infinity. In 

addition, the second term in the Kp integral does not con- 
! eur) 4 

tribute. Since a - 72 on the surface of the small sphere, 


the Kp integral yields 4amp. Hence, 


or=| 2 a¥e, 
Tpp' 


which is Coulomb’s law. 

If surface charges are present, the proof proceeds in a 
similar manner. However, the proof is not so simple for 
the case of point charges. 


7. DIELECTRICS, PHENOMENOLOGICAL TREATMENT 


We consider a charged condenser in vacuum formed by 
two parallel metal plates. We will assume that the areas 
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of the plates are large compared with their separation so 
that boundary effects need not be considered. Let the 
surface charge density be + on one plate and —w on 


hil 


+w —W 
te 


Figure 7.1 


the other (see Fig. 7.1). Under these circumstances a uni- 
form electric field exists between the plates, 


E,= +420, 


while the field outside is zero. The potential difference 
between the two plates is thus 


Pi— P2 = 4rood , 


where d is the plate separation. 

If one brings an insulator between the plates, then the 
following takes place: 

1. If the plates are connected toa battery so that their poten- 
tial remains constant, a current flows. Their charge in- 
creases. 

2. If, on the other hand, the plates are insulated, then 
their charge remains constant while their potential dif- 
ference decreases by a factor e: 


4mm 


Cit ae o- 


The factor ¢ is called the dielectric constant of the insulator 
and is, in general, greater than unity [A-2]. 
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The above phenomenon is interpreted as follows: A layer 
of charge of magnitude —(1—1/e) on the side near plate 1 
and +(1—1/e) on the side near plate 2 is formed on the 
outer surfaces of the dielectric. In conjunction with the 
charge on the plates, this results in a net charge of w/e at 1 
and —q/e at 2. 

The charge formed on the surfaces of the dielectric can- 
not be directly observed. Its existence can be demonstrated 
only through the potential difference. It cannot be altered 
by means of electrical conduction. Instead, it adheres to 
the dielectric substance and simply adjusts itself according 
to the external field. This type of charge is called polariza- 
tion charge to differentiate it from the true (conduction) 
charge which exists on the metal plates. The latter charge 
can be altered at will. 

In what follows, we will use 9, and w, to denote polariza- 
tion charge densities, o, and w, for true charge densities. 
Thus o=e,+¢@, and w=@,+w, represent total charge den- 
sities. 

If one imagines a thin region of vacuum to exist between 


Vacuum Dielectric 


Figure 7.2 


the condenser plates and the dielectric, then the relation- 
ships between the charge densities and the fields within 
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the condenser are (see Fig. 7.2): 


fp a Oe 


EY— EY = 4n0,, where w,=— (1 — ;) w. 
€ 
Hence 
: al 4 
EY = 4n|}o—(1—= = w, 
Fs 
or 
wk 
EY — - EY, 
€ 
and 


4nw, = — (e—1)E®. 


Figure 7.3 


the vacuum, in the absence of a true charge density, 


ele ek” =0, 
BO + EY = 4a0,=—(e—1)ED. 


This relation can easily be generalized to apply to the 


ny, Ne 


Figure 7.4 
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bounding surface between two arbitrary dielectrics by 
imagining a thin vacuum region to exist between them 
(Fig. 7.4). In this case, in the absence of true charge, 


&,E£,.+ &, E,, = 0 * 
If there is true charge present, 
€, HL, + & En, = 420: . 


For vacuum, «é is equal to unity. 
If weintroduce a new vector, the displacement D, defined by 


D=cE, = Ta 


then the most general relation for the bounding surface 
between two insulators becomes 


Da, + Dn, = 400; . [7.2] 
Since 
E,,+ E,, = 400 [7.3] 


(where w=w,+w,) is still valid, it follows that 
(e,: = 1)£,,+ (€.—1)E,, = — 4x0, . [7.4] 


This is an important property of matter. 

In the uniform field of the parallel plate condenser pola- 
rization charges are produced only at the surfaces of the die- 
lectric. However, in the most general case, a volume density 
of polarization charge 9, is produced within an insulator. 

Since the w,’s are the surface sources of the field D, then, 
by means of a limiting process (a series of stacked sheets, 
each with true surface charge density w,), it can be seen 
that the @,’s, the volume densities of true charge, are the 
volume sources of D: 


div D = 40, . [7.5] 
Because 
div E = 470 [7.6] 
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(where o=0,+¢0,), it also follows that 
div {(e—1)E} =—47o,. - er] 
If the medium is homogeneous, that is, if « is independent 
of position, then 
(e—1)divE =— 4x, . [7.8] 
If we apply Gauss’s theorem to the vector («e—1)E in the 


usual manner (excluding, however, the surfaces with w, +0), 
then we see that 


Jeoavt rd 0. [7.9] 


Thus, the sum of all polarization charges is zero. 

If w,=0 on a boundary layer between two dielectric 
media, then D,, the normal component of D, is continuous, 
while D,, the parallel component, is, in general, not. 
This results from the fact that HE, is continuous, while 

=—gradg, with g continuous, is still valid just as before. 


8. ELECTRON THEORY INTERPRETATION OF DIELECTRIC PHENOMENA 


In electron theory the polarization charges, introduced 
so far on purely phenomenological grounds, are traced back 
to electric dipoles. 


The electric dipole 


We first consider two point charges +e and —e separated 
by a distance d. The electric dipole moment of this arrange- 


ie 


T po. 
pa. 


—e +e 


Figure 8.1 
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ment is defined by the vector 
p— ea. 


(In general, for an arbitrary number of charges whose sum 
is zero, one defines 


P = 2 ess, 


where the x; are the position vectors of the charges. Since 
de, =0, this vector is independent of the position of the 


origin of the coordinate system.) 
The potential due to the two charges is, according to 
Hig 3... 


Now, if r>>d (i.e., either r— oo, or else €—> co as d—0 80 
that ed remains finite), this arrangement of two charges 
is called a dipole. A Taylor expansion of the potential due 
to this dipole yields 


i 
Pp = £(d- gradg) ==\(p-erad. as [8.1] 
pg Tro 


(Here, grad, denotes the gradient taken with respect to 
the coordinates of Q.) The field produced by the dipole is 
thus 


= = 
E, = (p:gradg) , 3 


PQ 


[8.2] 


Written in component notation, 


Up, — Tor ea 
i (x Di a) ae = (i, k = 1, 2,3). 


We now consider a finite piece of an insulator. Let there 
be a dipole present at each atom or molecule. (The ques- 
tion of how these dipoles originate will, for the present, be 
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avoided.) The field produced by these dipoles is then 


E, = > (p-gradg) ~~ *°. - [8.3] 
@ Te 

We wish to transform this sum into an integral by making 
the approximation of replacing the discontinuous distribu- 
tion of dipoles by a continuous one. Let p be the electric 
dipole moment per molecule and N the number of mole- 
cules per cubic centimeter. Then, the dipole moment per 
cubic centimeter, P= Np, is called the polarization. The 
field then becomes 


Ep= -| (P+ gradg) ~~, ~° dV. [8.4] 
Tre 

We wish to transform this integral. For our purposes we 

will need a generalization of Gauss’s theorem: For arbitrary 

vectors 


(A-grad)B + (div A)B = div (B, A), [8.5] 


where the divergence on the right-hand side is contracted 
with A. In component notation 


¥(4ean) Bet (SE) Be= Se (B,A,). [8.6] 


t 


For regular vector functions then, according to Gauss’s 
theorem, 


[(4-eraaypav— —[taivayB av -pA.B df. [8.7] 


In accordance with this relation, we obtain 


mae --[aiv TEN cad av — p Pome *8 te 
Tre Tre 

Here, P‘* represents the inwardly-directed normal compo- 

nent of P. This form for E, shows that our continuous distri- 

bution of polarization is completely equivalent to a volume and 
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surface distribution of polarization charge. Indeed, from the 
first integral 
0, =—divP, [8.8] 


while the second integral yields 
Wy = — PX” 


for the polarization surface charge density, provided that 
the region of integration lies only to one side of the surface 
so that the surface bounds the piece of insulator under con- 
sideration. On the other hand, the surface may be discon- 
tinuous and located within the insulator in such a manner 
that it is enclosed on both sides by the region of integration. 
The total polarization charge density on such a surface is 
given by 

Wy = — (Pa, + 4) . [8.9] 


This second case is the more general one and includes the 
first as a special case (P,, =0). 

Gauss’s theorem applied to the vector P shows once again 
that 


a Jeoav w,df = 0. [8.10] 


Furthermore, it is seen that in order to preserve the phe- 
nomenological laws of the previous section, P must be iden- 
tified with (e —1)E/4x. Therefore, 


D=E+4nP. [8.11] 


The dielectric constant ¢ can also be expressed in terms of 
atomic quantities by making the obvious assumption 


p=’, [8.12] 


in which « is the molecular polarizability and E’ is the 
effective field; i.e., the field prevailing at the position of 
the dipole p but not including the field produced by the 
dipole itself. If we make the rough approximation that 
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E'= E, then 
P= NzE 
and 
€=1+472N. [8.13] 

Erapirically, it is found that «> 0 for electrostatic fields. 
Comaequently, 2 ia greater than one [A-2]. 

We now with to refine this theory by considering the 
difference between E’ and E. The macroscopic field is just 
the average of the microscopic field: 


1 


hi we denote the field of a single dipole by e, then 


E= E'+ wfe-av. 


K 


Here, we have aseurned that the medium is isotropic. The 
integral extends over a sphere that contains the dipole. 
We will see that it is independent of the size of the sphere. 


y 


& 
Cylindrical coordinates o, 7 will be used. We represent the 


Calculation of |edV.* The geometry is shown in Fig. 8.2. 


P 


Figure %.2 


* LuwpBLaD, Ann. Phystk 57, 133 (1918). 
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radius of the sphere by R and P(e, z) is a point on its 
surface. Clearly, 
R?=0?+2?, 


From symmetry considerations, only the « component of 
La a 
= — A 
e e-gra (= ~) 


can contribute. Since 


oo 1 1 | 
e\Vert(C—aye Vert (C+ay’ 


the integral becomes 


é.=— 


RR +e=+V R92 


6) 1 1 
e.aV = —e[2nodo[ at; 7 —— See 
i Cet C—aye Ver+(Ctay 


0 —2=—V Ro? 


R 


~ | tao ae : —— 
: ie: (2— a)? es 


R 


=~ [treads z 4 
J vere Qaz VR+ a? rl d 


+R 


Anz dz 


——, as] ——— —___._ ——. , 
IVR + at 2az 


With the substitution w2=— R?+-a?— Yaz one obtains 


Rta 
2 — qp2 
AOU == = Oo thay Vee @, 
2a? 
x |R-a| 
tee : (ke? tL a?)u —443 R+a 
2a? |r—al e ’ 


If we now assume that R>a, then this expression becomes, 
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in fact, independent of R: 


4 4 
feav=—Zera=— 2p 
x 
Here, p= 2ea is the moment of the dipole. Thus, as the 


contribution of the field of a dipole to the macroscopic 
field, we obtain 


4a 


Kx 


It is essential that the above calculation be performed with 
a finite dipole (a-+ 0), since otherwise the integrand would 
diverge when expanded in powers of a. 

For isotropic media there results, then, 


E—E'——- Np=8—— FP. [8.14] 


In general, this formula is not valid for anisotropic media. 
As a consequence of 


al: 


and 


Pe Nake Na (e+ FP) 
alt 
= Nat +5le-n]E, 


the relation between the molecular polarizability and the 
dielectric constant is given by 


ma= jl ah 
wor eee! : 8.15 
e+2 Tica [ 


This is the Clausius-Mossotti relation for isotropic media. 
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If e—1<1, then e423 and 


e—lw4nNa, 


and we thus obtain the approximate formula of Eq. [8.13]. 


9. THE POTENTIAL PROBLEM 


On the basis of Coulomb’s law, we are in a position to cal- 
culate the potential due to a given charge distribution. 
However, in practice, other quantities are usually given and 
the charge distribution is ordinarily only implicitly deter- 
mined. 

The conditions which » must satisfy are as follows: 


1. At the boundary surfaces of conductors, either 


(a) —$e(Op/dn)df =4ae for an isolated conductor with the 
total charge specified. In addition, it is known that 
y =const. (The value of the constant is, however, not 
given in this case); or 


(6) p=, if the conductor is connected to a source of con- 
stant potential. In this case the potentialis given but 
not the charge. 


2. At the boundary surfaces of insulators, if w,=0, gm and 
edy/on are continuous. 


These are the boundary conditions for y. In addition, 
there is also the differential equation 


div (egradg) = 0, [9.1] 


which is valid for vacwum and for insulators where ¢ is con- 
tinuous. The region inside conductors is of no interest since 
we know that @ is constant there. In vacuum where ¢=1 
or, in general, where ¢=constant, this differential equation 
becomes 

Vig= 07 [9.2] 


It can be shown that under these conditions a solution 
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exists. The proof will not be given here—physically, the 
existence of a solution is obvious. In addition, it can 
easily be shown that the solution is unique provided that 
eé>0 (which is physically obvious [A-2]) and that gy at 
infinity vanishes at least as fast as 1/r. 


Proof of the uniqueness of the solution. Assumption: Two 
solutions g, and g, exist, each satisfying all the required 
conditions. We then consider the difference 


Ye fa 
At conducting surfaces, either 


) Fe(Ow/dn) df= 0, y= constant 
or 
(b) p=0, 


that is, $ye(oy/dn) df = 0 in both cases. 

At the boundary surfaces of insulators, eOy/On and yp are 
continuous, and div(egrady) =0 over the whole region 
where « is continuous. Because of 


div (ye grad py) = y div (e grad py)+ e(grad p)? , 


Gauss’s theorem, applied to the inside of a very large sphere 
from which all metallic bodies and discontinuous surfaces 
are excluded, yields 


reset pres l— hy (2 = + est) af 


conductor insulator 
a surfaces surfaces 
sphere 


=|elgraayy ay = 0, 


since all three integrals on the left vanish (the first because y 
vanishes as 1/r as roo). Because e>0, the integrand is 
positive definite. It can vanish only if grady=0 and 
y=const. However, since y=0 at infinity then y=0 
everywhere, thereby proving the uniqueness of the solution. 
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Mathematical remark. The differential equation 
div (egrady) = 0 


is equivalent to the variational problem 


1 minimum for given 
rp eae 2 — 0) == 
7 z, [ecerady) i deep: boundary conditions. 
The sum over é is to extend over all conductors on which 
the charge is specified (Case a). Thus, 


1 = i; festa grad (dg) dV — 26 dQ, 


iy ee 1 [ op a 
= ie div (grad) do dV + ig eae bea — Leeb.= 0. 


The surface integral vanishes on all surfaces except those 
of conductors of Case a. For such surfaces it yields deg, 
since n is the inward normal, and there remains ‘* 


1 
8J= — + | div (egradg) SpdV = 0. 


The differential. equation [9.1] is, therefore, a necessary 
condition. 

On the other hand, it can be shown that because «>0, 
then, together with the boundary conditions, it is also suf- 
ficient for a minimum. Indeed, let y satisfy the differen- 
tial equation and the boundary conditions. As the general 
solution of the minimal] problem we take 


P=9+ 9". 
We then obtain 


aft 
J(~) = J (py) + a e(gradg’)? dV, 


J(¢) > J (9), 


so that g is indeed a solution of the minimal problem. 
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Physical interpretation. The significance of the inte- 
grand which occurs in the variational problem is that of 
an energy density. Indeed, 


1 af i ut 1 
= 2). a) ee . — = Fe = : : 
gy o(arad ¢) a ee) ot ge) 


Here E?/8x is the energy density of the field in vacuum 
and (E-P)/2 is the potential energy per unit volume of the 
polarization. (The factor of 4 arises from the fact that the 
polarization is proportional to E.) The sum of these two 
thus gives the energy density in a dielectric. 


10. CURVILINEAR COORDINATES 


Expressed in orthogonal curvilinear coordinates u,v, and w, 
the square of the line element has the form 


3 
ds? = > (da,)? = ef du*+ e dv*+ edu’. [10.1] 


peal 


(There are no cross terms in the case of orthogonal coordi- 
nates.) Here, ¢,, €, and e, are functions of u, v, and w. 
The scalar product of a vector A with the line element is 


A,ds = A:dx = A,e,du+ A,e,dv+ Axe, dw. 


If one sets A=gradg, then the expressions 
oy Op oy 
— a —d 
cae i ee eis - 


and 
(grad q),€, du + (gradq),é, dv + (grad ¢),, és; dw 


must be identical; i.e., the wu, v, and w components of 
gradgy are 


10 10 1 op 
= 5 =i w=—z. [10.2 
(grad ~),= e, Cu’ (grad P)» es Ov b] (grad ~) es Ow [ 
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The surface element tensor is 
Ofun = €,0,dudv, df,,, = @e,dvdw, Afi. ¢,dudw , 


and the product of the surface element with the component 
of curl A perpendicular to it is 


curl, A df = (curl A),,¢e,e, dv dw + (curl A),e,e,; du dw +- ... 


Employing Stokes’s theorem, 


A dp (eke 2 cee ae 
ov ow 


one obtains, upon comparison, 


1 
€5 3 


O(ésAw)  O0(€,A,) 
ov ow 


(curl A), = 


| , ete. [10.3] 


It can easily be verified that curl(gradg) =0. 
Because 


A, df = A,e,e,dudw + A,e,e,dudw + A,,e,e,dudv, 


Gauss’s theorem can be written as 


re) 
p And ={[5 (e,¢€,A,,) i (e,e,A,) +. 


Since the integrand in the second term is div A dV, and 
since 


..|dudvdw. 


dV = ¢,¢,¢,dudvudw, 
then 


: ) 3) 
divA= las (€,€3.A jee ov (€,€,A ae Ow (€,€,.A .) [10. 4] 


€; €2 @s 
Upon substituting A= gradg in this last equation, the 
Laplacian operator V? = div grad is obtained. Thus, 


6) (a =| ) ie G5 4 4) ie op 


1 
Ou\ e, Ou ov €, Ov ar, eal [10-5] 


C1 Ose 


Vg 
” Ow \ €, Ow 
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This can also be derived from the variational principle: 


Jieraayy adv = min- 


11. EXAMPLES OF SOLUTIONS OF THE POTENTIAL PROBLEM 


By means of examples, we will illustrate two methods for 
the solution of the potential problem. 


a. Method of images 


1. A point charge e opposite an infinite conducting sheet. 
We assume that the sheet is grounded so that its poten- 
tial is gp=0. The method of images consists of assuming 
@ virtual charge —e located at the point which is the 
image of the position of e mirrored at the sheet. This 
image charge, together with the charge e, is to determine 


C0 


Figure 11.1 


the field to the right of the sheet but is to have no effect 
on the field to the left of the sheet (Fig. 11.1). Wethus have 


<0: op=0; «w#>0: dea 
This @ satisfies the equation V’m»=0 and is continuous for 
x=0. On the basis of our uniqueness theorem it is, there- 
fore, the solution to the problem. 
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To find the surface charge induced on the sheet we note 
that, 


é : é 
er (tp—a)—— (tp+ a). 


For «,=0, 
; Be 
re i, = — = = 400. 
Thus, 
—_ ea 
Ten 


The total charge on the conducting sheet is then 
€é =| w2r0do =| w2ardr=—e, 
0 a 
where r= V/02-+-a2, rdr=ode (0 is the distance from the 
origin to a point on the sheet). 


2. A point charge e opposite a semi-infinite homogeneous 
dielectric. In this case the boundary condition is 


ory 
= \ae). a 


where ¢ is the dielectric constant of the medium. We again 


Figure 11.2 


assume a mirror charge —e’ which contributes to the field 
to the right of the boundary, but this time the magnitude 
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of e’ remains to be determined. Then, for 


e 
e>0: =-——, 
, 


so that g is continuous at r=r’. Hence, 
op a oy a 
(32) =+S e+e, e(=2) = + Bo? )- 


The boundary condition thus requires that e+e’= e(e—e’). 
Thus, 
veo 


e= 
e+l1 


The solution of our problem is, accordingly, 


Me 5 


2>0: p=--——-— 


3. A point charge e opposite a conducting sphere. To every 
point P can be assigned a ‘‘mirror” point P’ by means of 
the relation (Fig. 11.3) 


Figure 11.3 
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Since triangles OP’Q and OQP are similar (they have a com- 
mon angle at O and o’/a=a/o), then 


r’ /r = ale = e'/a=V o'/e = constant 


for all points Q lying on the surface of the sphere. This 
fact can be used to solve the problem. 


Case 1. Grounded sphere. We imagine that a virtual 
charge —e’ at P contributes to the external field. Then, 


p= : -< outside the sphere, 


yp =0 inside the sphere. 


In order that » be continuous on the surface of the sphere, 
the relation 


must hold. From this, e’ is constant, and the problem is 
solved. The total charge induced on the sphere is —e’ 
(from Gauss’s law). 


Case 2. Insulated sphere. We specify that the sphere is 
uncharged. Thus, the image charge — e’ must be balanced, 
which we do by placing a charge +e’ at the center of the 
sphere. If & represents the distance from the field point 
to the center of the sphere, then 


e e e@' 

rs for ea, 
e' 

ae for Rae 


where e’ is still given by the formula for Case 1 above. 

For the case of the dielectric sphere, the solution cannot 
be given in closed form; instead, an infinite series is obtained. 
The same is true for two conducting spheres. 
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4. Conducting sphere in a uniform electric field. This problem 
can be considered as a limiting condition of Case 2 of 3 
above, if @ co and e->co in such a manner that e/o?—> E. 


Figure 11.4 


It is seen that a dipole of moment p is produced at the 
center of the sphere: 


p=e'o'. 
One can also begin with the following ad hoe hypothesis: 
x 
r>a@: gp=—Ha+A+B—, 
r<a: gp=A, 
where Ba/r? is the potential of a dipole. Defining «=rcos#, 
one can also write 
B 
r>a: g= 0088 (— Br+ =| +A. 
r 
The continuity of g can be achieved by setting the quantity 
in the parentheses equal to zero when r—a. Thus, 
B= Ka’. 


The constant A remains undetermined; it has, however, no 
physical significance. 
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The charge distribution induced on the sphere is 


4nw = — (=) = cost (2 + 4 = 3E cos?. 
Oni, : a? 


The total charge on the sphere is 


ovaj =o, since peosaf = 0. 


sphere 


5. A dielectric sphere embedded in a dielectric mediwm. Let 
the dielectric constant of the sphere be e; (=constant) 
and that of the surrounding medium be e, (=constant). 
At large distances from the sphere let the field be uniform 
and equal to HE. We make the hypothesis that for 


Cut % 
Tiree o.=—B(2 +z 7) =— Beoso(r—Z), 
r<a: 9; =—Dx=— Dreosd. 


The boundary conditions at r=a are 
x" 
ks B(1-4)=p, 


mer OMe, 
a = fae B14 wa) &= Dev. 


From this, it follows that the problem is solved if we take 


Semen - 3e, 
re 2&, 


It is interesting to note that for e;=« and «, =1, an 
expression similar to the Clausius-Mossotti relation of 
Kq. [8.15] is obtained, namely 


7! 
eer 


a’. 
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Because of this, it was possible for Clausius and Mossotti 
to explain dielectric phenomena by assuming that. a dielec- 
tric consists of small conducting spheres.® 


b. Application of function theory to two-dimensional problems 


if a spatial charge distribution depends only upon two 
of the cartesian coordinates, say # and y, so that the charge 
distribution is the same in all planes perpendicular to the 
third coordinate axis, then the field in all such planes must 
be identical and one speaks of a two-dimensional problem. 
In this case the problem can be solved with the aid of the 
theory of analytic functions. 

Let w be an analytic function of the complex variable 
e=a+ly: 

w(z)=u+w (u, v real). 

Then, because 


dw _ 


ou ov rl ov Ou 
Ox «Oy Ox oy 
From these it follows that 
o2 C2 
Vu=0 jand V%w=0, ~~ where Vi = Baa oye 


That is, uw and v satisfy Laplace’s equation. We will, there- 
fore, attempt to determine analytic functions so that their 
real parts u (or imaginary parts v) represent the desired 
potential. According to Liouville’s theorem, a function that 
is regular everywhere (including infinity) must be constant. 


’ See Enzyklonddie der mathematischen Wissenschaften, vol. V2, p. 329, or 
the original work of Clausius, Mechanische \W'drmetheorie (Braunschweig, 
1887-91), vol. 2, p. 62. 
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We will thus have to characterize functions according to 
their singularities. 

Since according to the Cauchy-Riemann equations, du/dn 
= 0v/ds, we have (Fig. 11.5) 


Ou 


ov : 
ae = bs, ds =[v] (change in v after one turn). 


If we identify the potential m with u, then 
xe ds = [0] = — Ar. 


Here, w is the charge in a cylinder of unit height. If we 
assume that the charge is distributed uniformly along 


a Figure 11.5 


a line of infinite extent normal to the x-y plane, then w 
represents the charge per unit length. 

Now, v is to increase after a complete circuit. This is 
a property of the logarithm of a complex number; thus 
w = In(z—Z)), v increases by 22. 

We obtain the solution for the case of a set of line charges 
perpendicular to the x-y plane with linear charge densities 
@, @,,... and intersecting the w-y plane at the points 
&y 2a, ..., if we set 


w = — } 2we log (2 — 2) . 
Q 
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The potential is the real part of w: 


Q=uUp= — > 2welogrpg, where rpg= Vige= XQ)? +(Yp— Yo)? . 
@ 


If the projection of the charge distribution on the z-y plane 
is a set of points, then the potential at these points has a 
logarithmic singularity. 


tee teeHeteHese + + 


Figure 11.6 


As a second example, we will calculate the potential due 
to a charged conducting wedge. Assertion: the potential is 
given by the imaginary part of the function 


Ws", 


where yw is real and must be suitably determined. For 
z=re® we have 


g=v0=Inw=rsinyd. 
Since the wedge is conducting, we must have 
u(r, P=0) = v(r', B= 2a — ax), all r, 7’, 
so that u(2x—a)=2. Hence, 


IU 


a) oe) 2 


The surface charge density w on the wedge is given by 


1 ov 
—4r0=-—| =pr'"cosw))| = pre. 
OT BB l ono ae 
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The total charge e(r) contained within a circle of radius r is 


2%—a 
| Pe Q— OH 


An 
— 4ne(r) =| 5,7 dd = — rh co8 u8| . = 12rt 
0 


12. MAGNETOSTATICS 


Electrodynamic phenomena cannot be described by the 
electric field alone; the magnetic field is also required. 
Just as electric fields interact with electric charges, magnetic 
fields interact with so-called magnetic dipoles, with whose 
aid magnetic fields can then be measured. (Later we will 
see that magnetic fields can also be measured with moving 
electric charges.) 

A very essential difference as compared with the electric 
case is that magnetic monopoles do not exist; they always 
occur in pairs so that the total ‘‘magnetic charge” is zero. 
An approximation to a monopole can be produced by 
removing one member of a pair to a sufficiently large dis- 
tance (i.e., a very long dipole). It was shown by Gilbert that 
Coulomb’s law is valid for such monopoles. 

In analogy to electrostatics we introduce the following 
vectors: 


magnetic field intensity H, corresponding to E; 
magnetization M, corresponding to P; 


magnetic induction B, corresponding to D. 


Aside from the very important difference that in magneto- 
statics there is nothing analogous to the true charge densities 
01, @,, the equations for the fields produced by magnetic 
dipoles correspond to those of electrostatics (Eqs. [2.4], 
[8.11] and [7.5)): 

cur = 0, [12.1] 
B=H+ 40M, [) 
div(B—0- [12.3] 
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At bounding surfaces, the components H, and B,, are con- 
tinuous. If we introduce a volume density o and surface 
density y of magnetic ‘polarization’ charge, then, cor- 
responding to Eqs. [8.8] and [8.9], 


o ——diyvyM=+_— divH, 
47e 


a 
4= — (Me, + Ma) = +7 (An, + He,) - 


In analogy to Kq. [7.9], it is also true that 


feaved rar=o. 


As in electrostatics (Eq. [7.1]), B and H are related by a 
constitutive equation, 


fo [12.4] 


where uw is the permeability and corresponds to ¢ [A-2]. For 
paramagnetic substances, w is greater than unity and is 
temperature dependent; for diamagnetic substances it is 
smaller than unity and temperature independent [A-2]. For 
ferromagnetic substances there is no unique relation be- 
tween B and H. (For example, B and M can be different 
from zero while H is zero.) 

Since H is irrotational, it can be represented by a scalar 
potential: 

H =— grady. 


Since B is source-free (solenoidal), it can be represented by 
a vector potential A: 4 


B=curlA. [12.5] 


However, A is still arbitrary to within the gradient of a 


4 Where there are no magnetizable substances present, the fields H and B 
are identical and A is then also the vector potential for H: H— curl A. 


56 ELECTROSTATICS AND MAGNETOSTATICS | Chap. 1 


sealar function. Thus, we will even be able to achieve the 
condition div A =0. 

As an example, we will determine the vector potential of 
the field of a dipole of magnetic moment m in the direction 
of the z axis. As we know from electrostatics, the field pro- 
duced is 

6) Bl Oo? 0? vel 
i m 5. (erad ;) = eas ” Oy dz’ a r 


Here, differentiation is performed with respect to the co- 
ordinates of the field point. 
Assertion: the vector potential for H is 


mx x aL 
Az = =— mx grad —. [12.6] 
Thus, 
af O01 
As = ees. eee ee = 
z Oy r? tf Ox 7? z 0, 
so that 
divA=0. 
For the components of H = ecurlA one obtains 
e 1 
T= = 
a ante r’ 
- o 1 
pee ee 
” By de yr’ 
o? oF. 1 oi 
H,=—m(4+5.\- =+me2 
‘ dat f = r ” Oar’ 


since V?(1/r)=0. Thus, we have verified that H can be 
represented by the A given above. Clearly, for a volume 
distribution of magnetization, 


Mx xp _ 
peer 


A= {aV, [aveMegraa, ~, [12.7] 
oP 


Yop 
where M, representing the magnetic moment per unit vol- 
ume, is, by definition, the magnetization. 
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It should now be mentioned that in addition to the 
previously employed obvious analogy between electric and 
magnetic quantities, 


H E 
B D 
M e 
be é 
there exists still another more fundamental one, namely, 
H D 
B E 
— M id 
1/u € 


as we Shall see later. 


13. UNITS AND DIMENSIONS 


The definition of units and dimensions is, to a certain 
degree, a matter of convention. The number of funda- 
mental dimensions can be reduced by arranging for certain 
quantities to be dimensionless. For example, either a new 
dimension can be chosen for electric charge, or else charge 
can be referred back to centimeters, grams, and seconds 
by choosing the constant k in Coulomb’s law, 


aa, 
72 


to be dimensionless: [k]=1. This definition is arbitrary. 
(In the analogous case of the gravitational law, 
M, Ms 

pe! 


Fa] 


f is not dimensionless if the dimensions of mass, length, 
and force are already fixed.) Setting [k]—1 corresponds 
to the cgs system. Thus, since V= ke?/r represents energy, 
the dimensions of charge become 


{[charge]— (erg-em)? = ({mass] [length]? [time]~*)* = (ml*t-*) 
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and 
[field] = [charge]em-? = (erg: em-*)? = (mI i)? 


The absolute electrostatic charge unit (esu) is defined so that 
k is dimensionless and equal to one when K and r are 
measured in egs units. 

On the other hand, if one sets k=1/47, then Heaviside 
units are obtained. In the following, these will be denoted 
by a subscript H. These units have certain advantages. 
It follows that 


Gm es field intensity = E, = i . eg = V4ne, 
where e is the charge in ordinary esu units, 

it 
V4 
Then the factor of 42 does not appear in certain other 
equations: 


E, E, K=<E =——_@hh 


divE,=o6,, instead of divE=4zo, 


while the energy is expressed as 


TC 


al 3 
W— — | H*avV=—— iz 
3 Hed 5 | Bea. 


Translator’s Note: Another common system is the rationalized mks or Giorgi 
system. This system has the advantage of the rationalization (elimination of 
the factor of 42 in many equations) offered by Heaviside units, and also has 
the consequence that potential, current, and resistance have the “practical’’ 
units of volts, amperes, and ohms. 

The fundamental quantities and units in this system are length (meters), 
mass (kilograms), time (seconds), and charge (coulombs). The coulomb is 
defined as the charge which, on each of two bodies separated by a distance 
of 1 meter, produces a mutual force of 1 kg:m/sec? (=1 newton =105 dynes). 

In the rationalized mks system k=1/47e,, where 


&) = 8.454 x 10-1? coul*- sec? -kg-!-m-> 


is called the permitivity of free space. Dielectrics are characterized by a per- 
mitivity «= ee where x. is known as the relative permitivity or dielectric 
constant of the material. In magnetostatics the permeability of free space is 
Hy = 42 x10-7 kg-m-coul-?. For magnetic materials, the permeability is written 
as #=xm,u, in analogy with the electric case. 


Chapter 2. Steady-State Currents 


14. THEORY OF STEADY-STATE CURRENTS 


a. Definitions 


The current J is the electric charge that passes through 
the total cross-sectional area of a conductor per unit time: 


_ de 
=a: [14.1] 

The current density iis defined as a vector in the direction 
of the current flow with magnitude equal to the amount 
of charge flowing through a unit cross-sectional area per- 
pendicular to i in one second. If the current density in 
a conductor is independent of position, then 


|i | — q ’ 
where q is the cross-sectional area of the conductor. Since 


the current through the surface element df in the direction 
of the normal n is given by 7, df, then, in the general case, 


r=fiaf. [14.2] 


Ordinarily, current is measured in electromagnetic units 


59 
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(emu). It can, however, also be defined in electrostatic 
units. In the latter case the dimensions are 


{J]= [charge ]sec = (erg-cm-sec~)? = (mi5t-)?, 
(i] = [J ]em-? = (erg-em-- sec)? = (mI)? . 


b. The continuity equation 


Because of the conservation law for electric charge, the 
amount of charge flowing out of a closed surface per unit 
time is equal to the decrease per unit time of the charge 
contained within the surface: 


de d 
} S eS ee - a 14. 


If the surface integral is transformed into a volume in- 
tegral in accordance with Gauss’s law, then, since the vol- 
ume under consideration is arbitrary, it is clear that 


) Saye 

~ ive 0%. [14.4] 
This relation is quite generally valid and is known as the 
equation of continuity. 


In the case of steady-state currents, all functions are, by 
definition, time independent: 6/ct=0. Thus, it follows that 


dive— 0, di. a); 07 [14.5] 
c. Ohm’s law 


Experimentally, it is found that when a current flows 
ina material, the following phenomenological relation exists 
between the vectors i and E: 


i=oE. [14.6] 
Here, o is the conductivity of the material in question. 
The above equation is called Ohm’s law. It can be ex- 
plained microscopically if it is assumed that a resistive 
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force proportional to velocity acts on the moving charged 
particles that constitute the current. 
Let N; and N_ be the number of particles per unit vol- 
ume having charges +e and —e, respectively. Then, 
oe = (Ni—N_)e, 
‘= N.:€°0,— iN Pee, 
where v; and v_ are the velocities of the positive and neg- 
ative particles. Now, Ohm’s law implies that each par- 
ticle is subject to a resistive force 
K,=—w,v,, 
where w is the specific resistance. The idea is that a sta- 
tionary state is established where 
K,=+ ¢E, 
me =—w,v,t+ceH=0. 
In practice, such a stationary state is formed in a very 
short time. Then, 


e 
= A= Ws. E, 
from which it follows that 
i=oE, where o=é (= +2) ; (aae7 | 
W+ w_ 


The first term in o is practically zero for metals. The fric- 
tional force on the electrons in a metal is proportional to 
their velocity with respect to the positive ion lattice. 
The dimensions of conductivity are? 
fol|=sec. 
Just as before, we have 
E =— gradg, curl E = 0, 
as long as the currents are time independent; that is, it is 
impossible to gain work from a charged particle which moves 
around a closed path. 


In electrostatic units. 
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For a homogeneous conductor (i.e., a conductor of homo- 
geneous material and of constant cross-sectional area q), 


Yi 2 


Figure 14.1 


Ohm’s law takes on the form familiar from experimental 
physics (see Fig. 14.1): 


V— 2 = L|E| 3 
J =qlil, 
so that 
i—-PHJsIk, [14.8] 


where R=I/(oq). The resistance Rk depends upon the phys- 
ical dimensions of the conductor. 

As we shall soon see, Eq. [14.8] is also valid for a con- 
ductor of arbitrary shape. In this case, however, R can- 
not be so simply expressed. 

The dimensions of resistance are ? 


[k]= [ocm]-!= em-!- sec—=1/velocity . 


d. Units 


Current: 1 electromagnetic unit (emu) = 3 x 10" electro- 
static units (esu)—10 amperes; 
1 ampere = 0.1 emu = 3 X10° esu. 


Potential: 1 esu=3 x10" emu = 300 volts; 
1 volt = 108 emu = 1/300 esu. 


Resistance: 1 ohm=1 volt/1 ampere = 10° emu 
=1/(9 x 107) esu. 


2 In electrostatic units. 


Sec. 14 | THEORY OF STEADY-STATE CURRENTS 63 


é. Joule’s law 


Heat is produced wherever electric current flows, and 
the heat developed per unit volume per unit time is 


VT E— ol. [14.9] 
The dimensions are correct since 
[oH?|= erg-(sec:cm’)-?. 

This law can also be obtained from microscopic consid- 
erations. The work done per unit time by the foree K= eE 
acting on a particle is 

K-v=e(v'E). 
Thus, the work per unit volume and time is 
e:(Niv,—N_v_) E=i-E. 


This work expended against resistive forces is completely 

converted into heat in the conductor. It is called Joule heat. 
The total heat produced per second in a homogeneous 

conductor of length / and cross-sectional area q is 


Q=i-Eq =J(y.—9.) =J77R. [14.10] 


We shall see that the expressions J(y,—9,) and J? are 
also correct for arbitrary conductors. 

The dimensions of this heat developed per unit time 
(= power= energy/time) are 


1 volt x1 ampere = 1 watt= 10’ erg: sec-!= 1 joule/sec . 


f. Current flow in conductors 
In conductors made of a homogeneous material, o is inde- 
pendent of position and we have 
=—grady, i=ok=—ogradg, 
divi=0, divE=0, V*p=0, 
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and thus 9 =0. In general, however, surface charges are 
present. Furthermore, 


at the surface. 
In conductors made of inhomogeneous materials one has 


div (c gradg) = 0 (o #9). [14.11] 


This equation is very similar to that for fluid flow. The 
lines of current flow are, therefore, essentially the stream- 
lines of hydrodynamic potential flow. 

Mathematically, the problem of current distribution is as 
follows. The potential difference between the electrodes, 


Figure 14.2 


Gi—2, is given. At the surfaces, o(dp/On)=0. In the in- 
terior we have the differential equation [14.11]. These 
requirements determine the current distribution [A-3]. 
Since o(dp/dn) is nonvanishing only at electrodes 1 and 2, 
it follows from Gauss’s law that 


Cpe. OM gn 
1 2 


For fixed geometry, if y is multiplied by a constant factor 
then the current also changes by the same factor. From 
this it is seen that for arbitrary conductors, 


1i—-G=IJIR. 
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We can now also prove Joule’s law in a more general man- 
ner. The total heat developed in a conductor is 


Q@ =|ck?dV=|i-Eav= ~[i-eraapay 


== +{paiviav—[aiv(gnav=—d pinay. 


Since 7, 40 only at the electrodes, we have 


QG= (Vi—- Pa) = J? R. 


If several electrodes are present a sum is obtained. 


15. THE MAGNETIC FIELDS OF STEADY-STATE CURRENTS 


For the present we can again think of the magnetic field 
as being measurable with magnetic dipoles. 

In his famous experiment Oersted found that the line 
integral of the magnetic field produced by a steady-state 
current, taken over a closed curve enclosing the currenf, is 
given by 


p H, ds = = EGS, [15.1] 


Here, c is the velocity of light, J the current measured in 
esu, and J,,,, the current measured in emu. For an infinitely 
long current-carrying cylinder, the radial and azimuthal 
components of the magnetic field outside the cylinder are, 


in polar coordinates, 


FA 
H.=0, H,=—<. [15.2] 


The relation between the sense of rotation and the direc- 
tion of current flow is shown in Fig. 15.1. This association 
is, however, necessary only if H is considered as a vector. 
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On the other hand, if H is considered as an antisymmetric 
tensor, then the association is not necessary. 


Current J flowing out 
of the plane of the page 


Figure 15.1 


As a consequence of its introduction in analogy to the 
electric field intensity, the dimensions of H are 


[| = |£] = charge-cm=. 
From 
[J] = charge-sec" 
it then follows that 
fo] =em-gee. 


Thus, ¢ has indeed the dimensions of velocity. 

We will consider the above relations as a fundamental law 
for steady-state currents. Corresponding te the equations 
for E, is 


Ba df = Ane curlE = 0, 
we now have the analogous relations for H: 
4 
pHas= Za, div H=0. [15.3] 


(The second equation for H requires that H ,=9, since 
6H, df =2xrH, =0.) Equation [15.1] is valid within a con- 
cylinder 


ductor as well as without: 


4 An [. 
pH.as= 24, = fia, [15.4] 
0 


F 
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where J, is the current flowing through the surface F enclosed 
by the arbitrary curve C. If the conductor is cylindrically 
symmetric with radius a, then we know that in the steady- 


Figure 15.2 


state the current density is uniform [A-3]. If the curve C 
is a circle of radius 7, then 


r\2 
J,= a(2) (GS a) ) 


Hie (r<a), 
Ht 4m 5 (ae ome 1 
2 ee Fi): 
- = (Z) 271 
H, = — - (r>a) 


Application: the solenoid. We will assume that the 
spacing between individual turns is small compared with 
both the length h and the diameter of the solenoid. For 
this limiting case we can roughly approximate the real 
field by one that is uniform inside the solenoid and vanishes 
outside. For the path of integration shown in Fig. 15.3, 


S| 
a 4 


Figure 15.3 


-.-— —. 


“ 


o=—— — .. aa 
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then, 


ae 
ficat =a pH.ds = Hh= "2 NI, 


where XN is the number of turns andd the current per turn. 
Thus, 


i 


<< [15.5] 


aly 


The magnetic field within the solenoid is directly propor- 
tional to the current. 

Just as to the equation ¢#,df=4ae there corresponds 
the equivalent requirement divE = 479, there exists a dif- 
ferential form corresponding to the integral relation $H,ds 
=4zJ,/c of Eq. [15.4]. According to Stokes’s theorem, 


4 
founntay = i, Af. 
Since this is valid for arbitrary surfaces, then 
47 : 
curl H = - b= 451 a [15.6] 


This is the differential equation which determines the 
magnetic field for a given steady-state current. 

Important Remark. We have seen that diveurla =0 for 
any arbitrary vector a. Thus, Eq. [15.6] is correct only 
as long as divi=0. This, however, is no longer true for 
the time-dependent case since, in general, dg/ot 40. Hence, 
the relation between i and H as expressed in the form of 
Eq. [15.6] cannot be correct in this case. 


Calculation of H with the help of the vector potential A 


The requirement divH=0 is automatically satisfied by 
setting 


H=curla. {15.7] 
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The vector potential A must now be determined. As already 
mentioned in Section 12, A is arbitrary to within the gradient 
of a scalar function. That is, for 


A'=A-4 grady 


with y arbitrary, H=curlA’ is also true. The function p 
can be determined by requiring, for example, that 


divA=0. [15.8] 


This is convenient for the case of steady-state currents. 
Because of the vector identity 


curl curla = grad diva— Va, [15.9] 


the vector potential must satisfy the condition 
: 47 
grad div A — V?4A = ae [15.10} 


(As employed here, the symbol V? represents the Laplacian 
operator expressed in cartesian coordinates and applied 
to each component of the vector a. In curvilinear coordi- 
nates the situation is somewhat more complicated.) Equa- 
tion [15.10] is satisfied if the following conditions are ful- 
filled: 


yp) ny = == 7 [15.11] 


Each component of the second part of Eq. [15.11] is anal- 
ogous to the expression V?9=—4zo of Eq. [3.5] which 
we know how to integrate: 9, = i| 0,4V,/r,9- Therefore, 


a J ley [15.12] 
c Tp@ 


We have still to show that divA=0. By twice applying 
the identity [4.10], 


div (fa) = fdiva+a-gradf, 
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we obtain 


since 1,=0 and divi=—0. 
Now, in order to calculate H from A we make use of the 
identity 
curl (fa) = f curla + gradfxa. [15.13] 
Then, 
1 i 
Hy =curl, A = < ferad, =x igdG, 


PQ 


since curl, does not operate on a Hence, 
alt 
ce} - 7 =a V6 


in which € represents the unit vector t= (Xp— Ai tag: 


Limiting case: a linear conductor 


For a conductor whose cross-sectional dimensions are 
very small compared with its length, the following mathe- 


Figure 15.4 
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matical idealization (which, however, is not always admis- 
sible) can be carried out. 

Let Q be the cross-sectional area of the conductor and ds 
the line element (see Fig. 15.4). Then, 


fiar=[iapias) =|} djas =[zas, 


(conductor) (conductor) & 


since i is parallel to ds. Then, 


, a [15.15] 


and 


_ od [Asx (xp—x%e) JF t 
H, = i F = _ [asexs : [15.16] 


Cc PQ 


This expression is familiar under the name of the Biot- 
Savart law. 

It is often said that the magnetic field produced by the 
current element Jds is 


t 
at 


“qH’’ — —dsx 
s 


J 
C 
However, the ‘‘dH”’ so defined does not satisfy the equations 
from which we started, since the current element does not 
satisfy the continuity equation. In reality there are no 
isolated current elements. Consequently, the Biot-Savart law 


rs 
if f 
Me 
0 r P(r,0,0) 


Figure 15.5 


12 : STEADY-STATE CURRENTS | Chap. 2 


can be applied only to closed current loops or to infinitely long 
conductors, the latter of which can be thought of as being 
closed at infinity. 

As a test, we apply the Biot-Savart law to an infinitely 
long linear conductor (Fig. 15.5). If the coordinates are 
chosen so that ds is parallel to the z axis, then 7, —a@, 
=r=constant and 


1 ag 
+o Zz— 
J dz J ie J1 
B= 7) SSS 
e } (+r) c¢ (r+2?)F\_. cr 


—O 


This is the formula which we have already obtained in 
Eq. [15.2]. 


16. THE EQUIVALENCE OF THE MAGNETIC FIELD DUE TO CLOSED 
CURRENT LOOPS WITH THAT DUE TO A DISTRIBUTION OF 
MAGNETIC DIPOLES 


We will first consider a transformation of the vector potential 
due to a closed linear current loop, which leads to the concept 
of magnetic moment. For this we require a special form 
of Stokes’s theorem. According to Eq. [15.13], for a vec- 
tor a which is independent of position, 


curl (fa) = gradfxa. 


Applying Stokes’s theorem (Hq. [2.5]) to the vector fa 
(where a is constant), we get 


pifa)-ds =[(graa jxa)-naf, 


a-pfds = a [nxgradyaf. 


Since a is an arbitrary vector, it must be true that 


pias =| nx gradf df [16.1] 
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For the vector potential 4,= (J/c)fds,/r 
relation results in 


J Hl 
A 7 [nx grad, — df. [16.2] 
c Tee 


this integral 


PQ? 


Here, only the bounding curve (i.e., the conducting loop) 
of the surface over which the integration is performed is 
determined; the surface itself is arbitrary. On the other 
hand, we have seen in Section 12 that the vector potential 
due to a dipole of moment m is 


at 
A, = mx grad,—. [16.3] 
Tre 
Thus, it is seen that the vector potential, and therefore the 
magnetic field due to a closed linear current loop, is equivalent 
to a surface distribution of dipoles with a magnetic moment per 
unit area of 
m=-n, 
ce 

where n is the normal to the surface. The choice of the 
surface is arbitrary, except that the bounding curve is deter- 
mined by the current loop. 


Special case: a plane current loop. As the surface, we choose 
the plane area enclosed by the loop. On this surface the 
distribution of magnetic moment is uniform. A further 
simplification results if the field point P is located very 
far from the current: rp, >a, where a is, for example, the 
diameter of the loop. Thus, for all points Q on the surface, 7p, 
can be replaced by 7r,,, where O can be any fixed point on 
the surface. Then, if F is the area of the enclosed surface, 


al 1 
A,=— whe es grad,-— = — Mx gradp—., 
C fen 


Teo 


Me Fae 
c 
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Thus, the magnetic field of a closed current loop is, at 
large distances, equivalent to that of a dipole of magnetic 
moment M as given by the above equation. 

What we have just done for linear current loops can also 
be done analogously for volume distributions of current, as 
we shall soon see. These equivalences were discovered by 
Ampére. It was also Ampére’s idea that in reality it is not 
magnetic dipoles but, rather, current loops which produce 
magnetic fields. In this theory the elementary quantity 
is not the magnetic dipole; instead, magnetism is attributed 
to microscopic currents (called Amperean currents) within 
molecules and atoms. 


We now consider a volume distribution of dipoles by a 
generalization of the preceding considerations. If the 
magnetic moment per unit volume (the magnetization) is 
denoted by M, then the vector potential of the field pro- 
duced is 


A=|Mx grad, —aV.. [16.4] 


We wish to transform this expression to 


A 


6 On ie 


Y 


aV+ Ake df. [16.5] 


This represents a vector potential produced by molecular 
volume currents i,, and molecular surface currents das Ma 
order to carry out this transformation we need an integral 
relation which follows from Gauss’s theorem. Because of 
the identity 


div (axb) = b curla—acurlb, [16.6] 


Gauss’s theorem of Eq. [2.9], with the help of a constant 
vector a, yields 


a-bx naj=¢p (axb)-naj=—a- feust bay. 
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Since a is arbitrary, 
pnxbaj= curlbdV. [16.7] 


With the aid of this integral expression and the fact that 
curl (fa) = feurla+ gradfxa, 


we obtain from Eq. [16.4] 


A =| Feu May —[eur (=) adv 


=| Feuimav— op 2*™ ay, 


r 


A volume distribution of magnetization is, therefore, equivalent 
to a volume current density i,, and a surface current density jm, 
where 

in = © CUTIE, [16.8] 


Jn =—OenxM. [16.9] 


The previously considered case of the linear current loop is 
contained in this result as a special case. (M constant in 
the interior; in passing to the two-dimensional case the sur- 
face current goes over into a linear current, and there is a 
magnetic moment per unit area rather than per unit volume.) 

The equivalence of the molecular currents i,, and j,, with 
a distribution of magnetic moment indicates the possibility 
of attributing magnetism to such molecular currents, as 
was done by Ampére. 

To differentiate it from these molecular currents, we shall 
designate the quantity previously referred to simply as 
“current” by the name conduction current and write the 
corresponding current density as i,. In addition, there is 
yet a third kind of current, the so-called convection current, 
which originates when a macroscopic charge is in motion. 
If o is the volume charge density of the moving body and v 
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its velocity, then the convection current density is given by 
i, =v. [16.10] 


It was shown by Rowland that the magnetic field produced 
by such a convection current is identical with that produced 
by an ordinary conduction current; that is, we have 


curl H = = (i, + i,). [16.11] 


The equivalence of molecular currents with a distribution 
of magnetic dipoles as discussed in these paragraphs now 
demonstrates the second analogy to dielectrics mentioned 
in Section 12. That is, the average value of the microscopic 
magnetic field intensity must be identified with the magnetic 
induction B, and not with H. This follows from the fact that 
for the steady-state case, 


(1) divB=0, 


(2) our B= A i = 4 i, iy) [16.12] 
Equation [16.12] follows from the relation B= H--1inM 
of [12.2], and from Eqs. [16.11] and [16.8] for curl H 
and curl M. 

The induction B has discontinuous tangential componegts 
if surface currents j,, are present. By applying Stokes’s law 
to the integration path shown in Fig. 16.1, it is found 


Figure 16.1 
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that 


. 4y 
BY _ Bo — = NX jn = 40M,. [16.13] 


Here n is a unit vector in the normal direction and M, 
is the tangential component of M at the surface. 

The vortices (curls) of the B field are thus given by 
all the currents present. Analogously, the sources (diver- 
gence) of E were determined by all the charges present, 
and E was the average of the microscopic electric field 
intensity. Therefore, in this deeper analogy B corresponds 
to E. However, B is a source-free vortex field while E (in 
the quasi-static case) is a vortex-free source field. (The 
operators curl and div are, therefore, also interchanged.) 

In analogy to D=E+42xP, the equation H= B—411M 
is, from this Amperean viewpoint, only an auxiliary rela- 
tion. At a bounding surface the normal component of B 
is continuous while H has a continuous tangential com- 
ponent. The more fundamental analogy between magnetic 
and electric quantities is, therefore, 


B E 
AH D 
—M | ae 


17. PONDEROMOTIVE FORCES 


Until now we have defined H as being the force acting 
on a unit magnetic pole. However, H can also be defined 
in terms of the force exerted on a current-carrying con- 


+ — 
Se 
Ve: Far away 


Figure 17.1 
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ductor element. The forces exerted by magnetic fields on 
conductors are called magnetic ponderomotive forces. 

The force acting on a unit pole at P (Fig. 17.1) due to 
the conductor L is 

K,= 7 [= Ss—*) dVo. [17.1] 
c Tre 

On the other hand, the magnetic field produced by the 
unit pole at a point Q on the conductor is 


According to the principle of action and reaction (action = 
reaction), the force on the conductor due to the unit pole 
at P must satisfy 


ko — kop 4), — — Ki, 


where K, represents the total force (i.e., integrated over 
all points Q). Thus, 


iG = 7 [ox HypdV, ° 


It now appears that the force exerted by the magnetic 
field H,, on the volume element dV, is given by (1/c) 
-(igx H,,)dV.. This consideration is, however, only heuristic. 
In the first place, the assumption of a unit magnetic pole 
is not very satisfying; nevertheless, it can be realized if 
the second pole is sufficiently far away. The chief reason, 
however, why the above consideration cannot be rigorous 
is that the principle of action and reaction is applicable 
only to the whole conductor, and the force on an individual 
volume element cannot be obtained from it. 

It turns out, however, that the force per unit volume is 
indeed given by 


1 
k=—ixH. [17.2] 
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This is a new fundamental law of nature. From this it fol- 
lows that the magnetic force on a linear conducting element ds is 


aK =2 dsx H. [17.3] 


From this, in turn, it can be shown that the torque experi- 

enced by a linear conducting loop carrying a current J and 

located in a uniform magnetic field H is the same as that 

acting on arigid magnetic dipole of moment M =(J/c) if ndf; 
F 


that is D=MxH. Ampére performed careful investiga- 
tions on the magnetic ponderomotive forces acting on con- 
ductors and thereby demonstrated the complete equiv- 
alence—also with respect to the effects of these forces—of 
closed current loops with dipole distributions. 

The above force law is also valid for convection currents 
i, = ov (where o is the volume charge density of the moving 
body and v its velocity). Thus, the total force K = | kdV 
acting on a body of charge e= i edV is 


K=“vxH. [17.4] 


This iscalled the Lorentz force after the Dutch physicist 
H. A. Lorentz. This relation shows that H can be defined and 
measured with the aid of moving charges, just as E can be 
defined and measured with stationary charges. 

The following force law is thus valid for the general case 
of an electric and a magnetic field: the force per cm? is 
composed of 


k, =oE (the electric force) 


Bid , [17.5] 
a — a ixH (the magnetic force). 


Here, iis the total current density, including convection cur- 
rents if present. This lawis quite general and is not restricted 
to steady-state currents. 
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The total force acting on a charge e consists of 


K, =cE 


and [17.6] 


Kegs de 
c 


The validity of this force law has been confirmed by the 
deflection ofcathode rays (electrons) by electric and mag- 
netic fields. It should also be remarked that the Lorentz 
force does no work since it is directed normal to the velocity. 
If H, is produced by a closed linear current loop (denoted 
by 2), then 
J, [ ds, x (x, — x) 
H,= Boe, [ Pree] 
2 
and the magnetic force exerted on a linear current ele- 
ment ds, of a second conductor (denoted by 1) is 


Ji J, { ds,x (ds, x (x, — x2) 
Ky = p e —. 


[17.8] 


Empirically, this is all that can be said (the current 
element ds, can be made movable). On the other hand, 
it is dangerous to omit the integral and to seek a fun- 
damental law expressing the effect of one current element 
upon another. Before the appearance of the Maxwell 
theory, this was indeed done from the “action at a distance”’ 
point of view. In doing this, however, one goes beyond 
what can be experimentally determined, since the contri- 
bution of the other conducting elements cannot be neglected. 
Accordingly, the result is not unique. The only thing 
that can be required of such a fundamental law is that 
upon integration, Eq. [17.8] must follow. For example, 
Grassman has found the “fundamental law’’ 


dK?” _ 


J, J, ds,x (ds, x Ny) x Xs 
= 3 » where ite 

Cc ft 
12 Tie 
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Equation [17.8] results upon integration over ds,. How- 
ever, it turns out that 


ake” 4 —dKe”, 


Nevertheless, the principle of action and reaction is valid 
for the whole conductor. There are still other ‘“funda- 
mental laws” which yield the same results. 

Historically, the oldest form is that of Ampére which is 
constructed on the basis of certain axioms: dK” is to 
have the direction of ny.=(x,—4%,)/r2, is to be bilinear 
in ds, and ds,, and is to be homogeneous of degree zero 
in 742, ds,, and ds,. From this it follows that it must have 
the form 


wag? ee 7 (ds,-ds,) + B(ds,° T12)(d8,° N,»)} : 


The coefficients « and 6 are determined from the requirement 
that 


p ak: ) en =dK?’, 


2 


in accordance with Kq. [17.8]. 

These ‘“‘fundamental laws’? completely lack physical sig 
nificance, since in the case of time-dependent currents 
it is not at all easy to say from which current element the 
field (which is, indeed, a measure of the force) originates. 
From the point of view of field theory, one cannot in- 
quire about such ‘‘fundamental laws.” 


18. THE PRINCIPLE OF ACTION AND REACTION FOR ELECTRIC AND 
MAGNETIC FORCES. THE MAXWELL STRESS TENSORS 


In this section, only volume distributions of current and 


charge are considered. 
We first introduce the concept of a tensor. By a tensor 
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in a 3-dimensional space is meant a system of 9 quantities 
T.,,(4, k=1, 2,3) which, under a transformation to another 
coordinate system, transform as the products a;b, formed 
from the components of two vectors a and b. (It is, however, 
not necessary that the 7,, be representable in the form a,b,.) 
A tensor can have special symmetry properties: if 


UEP Le the tensor is symmetric; 


Six =—Sriy the tensor is skew- or antisymmetric. 


An antisymmetric tensor in a 3-dimensional space transforms 
like a vector if only rotations (but not inversions) of the 
coordinate system are considered.’ Since S,, transforms 
as a:b,—a,b;, an antisymmetric tensor transforms as a 
vector product axb. 

Force is an ordinary (polar) vector, as is the current den- 
sity i. Thus, according to Eq. [17.2], His obviously an anti- 
symmetric tensor. This can also be seen from the Biot- 
Savart law of Eq. [15.16]. Actually, one should write 
H,;, Hz, and H,, for H,, H,, and H,, respectively (with 
H,, =—H,,). On the other hand, the electric field intensity 
is an ordinary vector. 

We denote the scalar product of the tensor 7 with the 
unit vector n by the vector 


I, == CTs nae de) 
in which 


Tui > Tim (¢ = 1, 2,3). 
k 
It is immediately seen that the 7, transform as the com- 


ponents of a vector. 
Gauss’s theorem has the form 


or; 
ib at =$ 3 Tamas, [18.1] 


* See the discussion on pp. 11 f. 
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which can also be written as 


fai have ) Td, [18.2] 
where 
Ode 


ox k 


, ro 
2 ; > a : [18.3] 
The principle of action and reaction is satisfied for electric 
as well as for magnetic forces if k, and k,, can be repre- 
sented by 
ivi, ke=div 1, [18.4] 
where 7 and 7 vanish sufficiently fast at infinity, so that 


Je aV= prea 0 


and [18.5] 


[roav— p T” df +0 


as the surface bounding the region of integration becomes 
infinite. 

We will now prove that the principle of action and reaction 
is satisfied for steady-state currents and stationary charges 
(i.e., for time-independent fields). For this we employ the 
previously derived relations of Eqs. [17.5], [3.4], [2.4] [15.6], 
and [15.3]: 


h, = oR = divE-E, curl E=0, 
4a 


i i! ‘ 
kn, = -—ixH= —curlHxi, divH=0. 
c 47 
(We use the equations corresponding to the vacuum since 
we are considering the microscopic fields.) These can also 
be written as 
i) 


= ve {E divE + curlEx E}, 


hy = (HdivH + curl Hx H}. 
“4 
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Consider the purely mathematical transformation 


(curl Ex E], = (curl, E)E,— (curl, E)E, 


dE, oH,\_, (dE, oF, 
“(as sat) Bs a ae 


= 


Ey ene 


kat 08> 


In the last expression two terms whose sum is zero have 
been added. For an arbitrary index 7 we have 


OX; Ox 


3 
[curl Ex E]; = » (Gn se) Be, 


from which it follows that 


= oH,\ , (OF; OF; 
a [B a Bp ( os = B,| 


6) é wt 


Here, H? is a scalar which has been obtained from the con- 
traction of a symmetric tensor. With the help of the 
Kronecker delta symbol 


Q” for t-2k, 
bn = [18.6] 


1 for i=k, 


which is itself a symmetric tensor with respect to orthogonal 
transformations, the tensor 


4aTy = = H;E,— dbx +E? 
can be formed. Then, it is indeed true that 


0 oe 01 2 
72 a, gg eB) 3a, 22": 


k 
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Equations [18.4] are thus satisfied for 


TS = _ {. B.—5 ba ; 

Wet (x, 1 sair| | [18.7] 
4 aoe 

These are the Maxwell stress tensors. If all currents and 

charges are restricted to a finite region of space, then the 

fields vanish at least as fast as 1/r? and the surface in- 

tegrals in Eq. [18.5] of 


i. 


@ 1 ey 1 2 


and 


To z |HUH-n) — nH] 
approach zero as the surfaces become infinite. 

The Maxwell stress tensors are symmetric. This implies 
that the resulting torque is zero (conservation of angular 
momentum). Proof: The torque per unit volume isd =xxk. 
Then, 


ten OD 
dj; = —4;,;= 0k; —2,k,= E(«. Oa, j ma 


& 
0 
aoal 2 a5, (@: 2 j,—%;T x)—Ty4+ Ts. 
k 


With 
Ly Ly= d:,, 


Gauss’s theorem yields 


faav=paxr,aj+{ av. 


The first term vanishes, and if 7,;=7;:, so does the 


second. 
A scalar can be obtained from any tensor by taking the 
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trace: 
T =) 7,7 scalar. 


(Since 7, transforms like a,b,, then T transforms like 
> a:b; =a-b, which is a scalar.) In the case of the Maxwell 


z 


stress tensors, 


t 


é ul 3 1 
S P= (P—5 P| --Z8=— ey [18.8] 


2 70 


al 


LIP=-—=—We, | [189] 


where W, is the energy density of the electric field. We 
shall see later that W,, is the energy density of the magnetic 
field. 

The Heaviside units are quite analogous in the electric 
and magnetic cases: 


1 
Cn — V 4x6 E =——E 
A 3 a ae ’ 
k, = 0E = OnE; , div E, = Ox; 
— 1 
a b] H \/ de b 
JL ita gs 
(P= 3x = 5 ieX He, curl H enn 
W.= }E, eee 


Translator’s Note: In the rationalized mks system (see note, p. 58), the 
equations differ by constant factors from those given above. For example, 
Oersted’s law (Eq. {15.1]) reads 


fH -ds=J, 
while the Biot-Savart law (Ha. [15.16]) is 


L [tox er, 
Hp= res a Ue where Tpg=Xp— XQ. 
PQ 


The current is measured in amperes and the units of H are ampere-turns/m. 
The fields B and H are related by a constitutive equation which can be 
written in several ways: 


B= uH= xmpH =p (1+ xm)H= »H+M, 
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Here, u.=42* 10-7 kg-m-coul-* is the permeability of free space, uw is the 
permeability of the medium (in kg-m:coul~*), xm- “/u igs the relative per- 
meability (dimensionless), xm=*mn—1 is the magnetic susceptibility (dimen- 
sionless), and Mis the magnetization. The units of both B and M are webers/m?, 
where 1 weber = 1 kg-m*/coul:sec. 

In the rationalized mks system, the following equations hold: 


Electric Magnetic 
ke= 0E, kn =iXH, 
divE=o, curl H=i, 


7 «= 48E*, Wmn= 4H. 


Chapter 3. Quasi-Static Fields 


The fields to be considered in this chapter are assumed to 
change but little during the time required for light to 
traverse a distance equal to the maximum dimension of 
the body under consideration. Thus, the finite velocity 
with which the fields propagate need not be considered. 


19. FARADAY’S LAW OF INDUCTION 


The essential law which, in the case of quasi-static fields, 
must be added to the previous equations is Faraday’s 


law of induction. When there are no magnetizable bodies 
present, 


ds 


Figure 19.1 
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the most general form of this law is 
id 
E,ds = — —- — : ; 
$ ,ds 5 5 | Heal [19.1] 
o F 


In contrast to the static and steady-state cases, the line 
integral of the electric field intensity over a closed path no 
longer vanishes. The minus sign in Eq. [19.1] can be 
interpreted in terms of energy: switching a current on re- 
quires the expenditure of work (see Sec. 20). 

Equation [19.1] can be transformed into 


1d 1 (oH, 
fount. bap=— 25 [H.ay=—2 [% df. 
F 


F F 


Since this holds for arbitrary surfaces, it must be true that 


P| 5 on 


— [19.2] 


This is the differential form of Faraday’s law of induction. 
For magnetizable bodies, B appears in the place of H: 


1d 
pe. ds = — aul B df, [19.3] 
c P 
10B 
curl E = — mar [19.4] 


From Eq. [19.4] it follows that 


divB —0. 


2] 


This agrees with our stronger postulate of Eq. [12.3], 
divB=0, [19.5] 


or, for vacuum, 
dia 0; [19.6] 
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If, at any given time, Eq. [19.5] or [19.6] is satisfied, then 


it follows from Faraday’s law that it is satisfied for all times. 
Furthermore, Eqs. [3.4] and [16.12] are valid. Thus, 


divE = 4xo [19.7 
and 
eurl B = = rs [19.8] 
or, for vacuum, 
4 
curl H = = a [19.9] 


The last two equations are approximations since from them 
it follows that divi=0, which condition, however, is not 
fulfilled exactly. All other equations in this section are 
exact. 


20. THE ENERGY OF CURRENT SYSTEMS 


We are now in a position to calculate the energy of a 
current system. The energy is equal to the work done by the 
induced current after the applied voltage is switched off. 
Since the force per unit volume is given by Eqs. [17.5], the 
work done by the system per unit time and volume is 


A= ok [20.1] 


We assume that the external field is switched off at t=0. 
Then, the work subsequently done by the induced field is, 
with the use of Eq. [19.9], 


w=fat{aar ~aifi-Eav— j, {ae Bounty. 
0 0 


Because of 


E-curl H— H-curlE =— div ExH 
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and Kq. [19.2], we have 


ioe) 


ye — ial H-curlE dV = if H-Hdav 


t=<o 
= — — af Hs i av ‘ 
leno 
(The surface integral can, in the quasi-static case, be ne- 


glected since the fields at infinity vanish sufficiently fast.) 
Thus, since H=0 at t= oo, we obtain 


1 
—_ — 2 
Ww a av. [20.2] 


This is the energy of a current system. This is also equal 
to the work which must be done in order to switch on 
the current. The energy W is positive because of the 
minus sign in Faraday’s law of Eq. [19.1]. 

Since divH =0, it follows that 


H=curlA [20.3] 


quite rigorously. On the contrary, it follows from Eq. [19.9] 
that 
divi=0, [20.4] 


and, as shown in Section 15, that 


ye | ay [20.5] 
cl; 
and 
divA=0. [20.6] 


These last three equations are only approximations for 
slowly varying fields. The energy W given by Eq. [20.2] 
can be transformed as follows: 


= a {Har gq | Hem aay. 
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After partial integration, dropping the surface integral, and 
using Eq. [19.9], 


w= a A-curlH dV = — fA ‘tdV. 


Thus, with Eq. [20.5], 


Up ip: 
= sa {|r aveave. [20.7] 


Equation [20.7] exhibits a definite analogy to the energy 
of the electric field given in Eq. [4.5]. The expression of 
Eq. [20.7] can be written for individual circuits as follows. 


a. One circuit. Here, 
Wass [20 8] 


The quantity Z, called the coefficient of self-inductance, 
is defined so that 


iy: se 
= = { [= aV,aV,. [20.9] 
ce Tpp' 


This is reasonable, since L depends only upon the geometry 
of the circuit. That is, if the current J is multiplied by 
some factor, then i at every point is also multiplied by the 
same factor. Here, it is only presupposed that the de cur- 
rent distribution can be used in the integrand. For rapidly 
varying fields, LZ becomes frequency dependent. 


b. Two circuits. In this case the energy of the magnetic 
field takes the positive definite quadratic form 


W= 31,014 LydiJ.+} Ld}. [20.10] 
Here, 


il incip 
Lyd, J, = nave fave nae [20.11] 
c Tpp' 
1 2 


If the conductors are filamentary, then 


iiidaed 
Lae I Ie PP pee? [20.12] 


Nie 


pees 
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This transition is allowed since r,,, remains finite. The coef- 
ficient of mutual inductance then becomes 


i ds,-ds, 
Ly, = L, = a pp =: [20.13] 
1 2 


On the other hand, the transition from spatial to filamentary 
conductors cannot be made for the case of self-inductance; 
otherwise, completely erroneous conclusions may result. 

The generalization to more than two circuits is obvious. 


Coefficient of self-inductance for a solenoid. According to 
Eq. [15.5], 
4nN J 
Usha oars 
Thus, 


1 1 1 1672? J? 
= ane 2 ial 2 = 
5 Ld =F dV =eerath 


1 4a%a?N? 1 (2naN)? 1 A? 


ace. 6k Tak? oe 


where A = 2zaN is the total length of the winding (see 
Fig. 20.1). 


batat 


Figure 20.1 


Dimensions of self-inductance. In cgs units, 


5 Shei a el Poe 


——i NAOT 6 


eau ? 
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In the volt-ampere (rationalized mks) system, 


volt - see 
1 henry =) _.__ _ — 10 em. 
amp 
21. TIME-DEPENDENT CURRENT FLOW IN CIRCUITS 


From Eqs. [19.2] and [20.3] it follows that 


10A 
curl (E+ Ae =(). 


We thus have the following generalization of the electro- 
static relation of Eq. [2.3]: 
E=— ~~, — rad 9. [21.1] 


We now wish to formulate an important relation for the 
variation of the current in a circuit. Let the current distri- 
bution in space be given. We assume that it is the same 


7 


Figure 21.1 


as for a steady-state current. We first consider an ordinary 
conductor with two electrodes. As demonstrated in Section 14, 


i=ocE, [21.2] 


while the heat produced per unit time is, on the one hand, 
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given by 
= y2 
Q -{= dV —J:R [21.3] 


from the definition of the resistance R, and, on the other 
hand, by 


Q= i-Eav = —[i-gradpav— sid dV 
because of Eq. [21.1]. If we neglect diva in 
div (yi) = p divi-+ i-gradg 
and substitute A from Eq. [20.5], we get 


,_ die 
= a. 1 ? ot 
=~ vinaf— aff T pp’ dV,daV, 


Differentiating Eq. [20.9] with respect to time results in 


a a 
ey — 2— 1 ft aV,dV, 


pp’ 


(since interchanging P and P’ does not affect the integral). 
Thus, we obtain 


Q=J*R=— pind) — LIF, Rai ed , 
where 
E=g'— 
is the electromotive force. We thus obtain the important 
equation 
LJ +RJ=e. [21.4] 
This could have been obtained more easily by assuming a 


filamentary conductor. The latter procedure, however, is 
not correct. 


96 QUASI-STATIC FIELDS | Chap. 3 


For a circuit containing a capacitor, an additional term is 
obtained: 


— poi. af = ed — (gy — G2") - 


The sum of the charges on 1’ and 2’ (Fig. 21.2) remains 
constant and equal to zero. What flows in at the left 


A +€ —e 


®, 


& 


Gr Px 


Figure 21.2 


flows out at the right; it is as if the current flowed 
without interruption. If the charge is denoted by e and the 
capacitance by C, then 


Vim Pa 


Qle 


On the other hand, 
Ce 
In this case the differential equation becomes 
Liao woe 5 = [21.5] 
or, differentiated, 
LJ Wa =é, [21.6] 


For the case of two cirewits, Q is calculated separately for 
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each. However, A is now due to both circuits. Thus, 


o -" ea fi dav, 


1 


sri iee dV py aie > 


Analogous to the preceding calculation one obtains 


Lip eines +EL+ GHA, 
o [207] 


. . € 
Ded, + Lyd, + Red2+ a = &2. 
2 


The terms e;/C; drop out if there are no capacitors in the 
circuit. The energy is 


1 1 if i 
B= sh Sit luAhat+sbatsa tag: [21.8] 


Because of 
J,= 41, J,= és, 
it follows that 


Hah t oh BA BS, [21.9] 
which is equal to the work done by external forces minus 
the heat produced per unit time. In the above differential 
equations, R has the character of a frictional force since 
the equations are not invariant when J is replaced by —J 
and t by —t. 

We have obtained ordinary differential equations rather 
than partial differential equations because we have assumed 
that the current distribution is known, and also because 
radiative effects have been neglected. 
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Integration of the differential equations for some interesting 
cases 


1. One cirewit without capacitance (an RL circuit). The 
differential equation is 


LJ +RJ=e.. [21.10] 
Switching off the emf, 
Sea Or t-2 0, 


6=0" for t= 0, 
gives the solution 


10% tO: J=h=3 (Ohm’s law), 
for t>0: G==55, 
ee 
J =dye * 


(The constant of integration has been chosen so that the 
current is continuous at t=0.) This solution corresponds 
to a gradual decay of the current. 


Switching on the emf, 


implies 
for t<0 7 =, 
for t>0 +7 (¥-2) =o, 
g-~ Ba o.6%, 
C_ 


so that J is continuous at t=0: 


J=2(1-e 


, 
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Thus J—>e/R as too. Initially, the current increases 
gradually because of the self-inductance. The final current 
is established in a time of the order of magnitude of L/R 
(proportional to 1/2). 


2. Applied electromotive force with periodic time dependence. 

For linear expressions it is convenient to use complex expo- 

nential notation. The real (or imaginary) part can be taken 

later, as required. That is, if we find a complex solution of 

a linear differential equation with real coefficients, then 

both the real and the imaginary parts are also solutions. 
We therefore assume that 


— jt 
E=(9e" 5 


where w is the angular frequency, equal to 27 times the 
frequency and try 
Ae, 


For this type of solution the differential operations are 
equivalent to certain multiplications: 


—~iw, 


dt 


2 
<~ (io) = — a. 


Equation [21.10] then yields 
(lin + RJA = &, 


Eo E(k sal iw) 


a rr, Re wir 


If & =|e|e%, then we can also write 
A=(|A| eft , 
where 


|A|=|e)|/(R?-+o2L)! and tand =—oL/R. 
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The phase of the current is shifted by an amount # with 
respect to the phase of the applied electromotive force, 
since 

Re (€) = | € | cos (wt + a) 


and 
Re (J) =|A|cos(@t+a+#). 


3. One circuit with capacitance (an RLC circuit). The dif- 
ferential equation is 
“ = ol 
Id + Rd + gee’ [21a 
a. Forced oscillations: ¢=«ge'*'. We assume that e<=0 and 
that ¢ is real and positive and try 
J nem 

Thus, Eq. [21.11] yields 


(~ wD + iok +3) A=&,io, 


1WE 
AS é = lA lee 
PST EP EUs aml 
where 
FS 1 if 
A) Sapna tand =~ 5 ( al 
AI V R?+ (ol —1/Cw)?’ R\°" Co 


The maximum of [Aj occurs when the denominator is 
minimum: 
ul 


on Com My 
1 
OS SSS 
VTC 
resonance. 
| Allie = bY 
R 
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The quantity tan? changes sign when going through the 
resonance. 


|A| 


Ww @ 


m 


Figure 21.3 


b. Free oscillations: e=0. Put 
J = Ae, 


The oscillation frequency w, is now not established by the 
emf but, instead, is determined by the differential equa- 
tion of [21.11]: 

ik 1 ai? 


: 1 
— Loy + Rion + G=9, ss alg 5 nC a: 


There are two possibilities: ? 


Re 1 L 
= i 1G! n<al/5; 
a R? 
— a, — —_ 
Re (J) = Ae cos( Gi Ea 


1 Translator’s Note: In addition to the damped and overdamped cases 
discussed in « and §, there also exists a third possibility, namely, when 
R?/4L?=1/LC. In this case, 

@) = iR/2L =i/+/LC 
and - 
Re(J) = Ae~(BI2Dt = get VEO _ gent , 


This is called the critically damped case. 
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(with appropriate choice of the time origin). This corre- 
sponds to damped oscillations. 


(B) eee /E. 
i ee R>2| ; 
ee) oe 

=N\or =) an Le)? 


Rel) =A-e (ee? ae 
Here, there are two solutions, corresponding to both signs 
of the square root. This is the aperiodic (overdamped) case. 
In the periodic case «, the oscillation frequency wp, is 
somewhat smaller than the resonance frequency @,n. 
If R<2V/L/C, then 


Wy—> Wm =1]// LC , 
which is sometimes called the Thompson equation. 


4. Two RLC circuits in parallel. If the mutual interaction 
between circuits 1 and 2 is neglected, Eq. [21.7] differ- 


Ts; i 
Lh, ies 
Figure 21.4 


entiated with respect to time yields 


. oe. 2 J ms 7 dle 
€=L,d,+ Ridit+ G = laht+ Bedi t G- 
‘1 2 
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In the spirit of forced oscillations we set 


J, = A, ert ’ 
J, = A, eimt ; 
Een 


Then, 

; ' 1 ; al 

tWéy = A, (- Ly? + Rytiw + Z| = A, (- Da? + Reto + z) ; 
1 2 


i 1 
0 EE a ee ee 
it Ay ive + R,iw+1/C, a ee 
J = iwe | zt [. i 
i —L,0?+Ryiw+1/0, —L,w?+R,iw+1/¢,}° 
For the case of direct current, in the absence of capacitors, 
1 if 
ie (z + x) . 


22. THE SKIN EFFECT 


Until now, we have assumed that the current distributions 
are the same as for the direct current case. However, with 
increasing frequency the self-inductance tends to restrict 
the current to the surfaces of the conductor. The resistance 
is thereby increased since the effective cross-sectional area 
is diminished. 

We will perform the calculation for a semi-infinite con- 
ductor. (Bessel functions appear for cylindrical conductors.) 


¢ Conductor a 


Figure 22.1 
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Let the current distribution be parallel to the y axis. (For 
direct current, it would be uniform.) We employ the equa- 
tions 


curl H= = ie [22.1] 
Ris 

curl E = — oH: [22.2] 

i= ee. [22.3] 

divH=0. [22.4] 


Because of Eqs. [22.3] and [20.4], it is also true that 
divE = 0. [22.5] 
From the first three equations it follows that 


4x Oi - 4n OE 


curl curlE = — ao a ae 
4n OH 
1 =——o—. 
eurl curl 0 ao Oy 


Because of the identity given in Eq. [15.9], consideration 
of [22.5] and [19.6] leads to 


4n OE 
meme [22.6] 
4n OH : 
VES dens: . 


These two relations are very similar to the equation for heat 
conduction. They are not exact because the usual ap- 
proximations for quasi-static fields have been made. Nev- 
ertheless, they suffice to explain the skin effect. We will 
assume that the fields have the form E£,(z), H,(z), all other 
components being zero. This hypothesis is sufficient. Let 
the fields be periodic in time: 


Ey = EX(eye™ . 
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Substituting into [22.6], 


@E 4now. 
a c se 


With k? =47owi/c?, the solution is 
K=Ac. 
Introducing the quantity x=+/2z0w/c, so that 
M21, k=+x(1+2), [22.7] 
the solution becomes 


H, = J e%% ei(xetot) 
and 
Re H, = Ae* cos (wt + xz). [22.8] 

We must restrict ourselves to the positive value of k since EB, 
cannot increase with increasing depth. Thus, E, decreases 
exponentially as z>— oo. (For w=k=0, one obtains at 
most a linear dependence of E, on z.) 

The effective depth to which the alternating current pene- 
trates, 1/x, is 


1 il 
% 27CemaW : 
since 
o/0*= Oca 


Numerical estimate. For a frequency v= w/27= 50 sec", 
for copper where o,_, =5.9:10 sec:cm-?, 


1jz~1cm. 


For the resistance of a cylindrical wire, a calculation yields 
the asymptotic formula (xa >1) 


ua 
R=kh>- 


Here, R, is the de resistance and a the radius of the wire. 
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23. THE LAW OF INDUCTION FOR MOVING CONDUCTORS 


The following considerations represent a preliminary step 
toward the theory of relativity. 

Let the conductor move with velocity v. The force on a 
moving charge is not eE but, according to Eq. [17.6], 


(E+ 5 ved). 


For a moving conductor, Ohm’s law must be modified to 
the extent that 


E* = E+ ~ (vx H) | [23.1] 


appears in place of E, so that the current density is given by 
P= ol. [23.2] 


(The above is correct only to first order in v/e. An exact 
consideration shows that o is dependent upon direction.) 

Now, purely mathematically, without any additional 
physical assumptions, it can be shown that if 


1d 
pb.ds=— 25 [Hal [23.3] 


is valid for a stationary surface, then 


*ae—-—ia 
Hts = a [aa [23.4] 


> 


holds for a moving surface. (The symbol — is to denote that 
the surface and its bounding curve are in motion.) 

This, together with the physical law of Eq. [23.2], indicates 
that for the same relative motion between conductor and magnet, 
the same current is induced, whether (a) the magnet is moved 
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slowly with respect to the stationary conductor (Fig. 23.1a), 
or (b) the conductor is moved with the same—but oppositely 
directed—velocity with respect to the stationary magnet 
(Fig. 23.1b), since }H*ds is the same in both cases. This 


Stationary 


— 


Stationary ~~ > 


Figure 23.1 


fact is more extensively generalized in the theory of rela- 
tivity. 

Proof of this mathematical theorem. For an arbitrary vec- 
tor B at the time ¢, it is true that 


d OB, d 
= [2.ar= ap df + = [Batovar. [23.5] 


vdt 


Figure 23.2 


We now construct a closed surface by supplementing the 
surface at times t and ¢+dt with the area of the cylinder 
swept out by the conductor during the time dt. On this 
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cylindrical surface, 
df = dsxvdt. 


Applying Gauss’s theorem to the closed surface gives 


fo pav—/2. oat |B. oar +[ -(dsxv dt), 


t+dét 


where the argument of B is evaluated at a fixed time t. 
Transforming the last integral by means of Stokes’s theorem 
and using dV=d/v,dt results in 


fasmar—[2.may— div BdV — at/as: (vx B) 
t+dt t (3) 


= dt}{v div B— curl (vx B)}, df. 
Substituting this into Eq. [23.5] yields the mathematical 
identity 


d 


ai B,, a=|{S + vdiv B— curl (wxB)} df. [23.6] 


This relation is valid for any arbitrary differentiable vector B. 
We now apply the above to the case B=H. It is es- 
sential to use the fact that 


divH=0. 


There then remains 


14 facay= 2] (28) ap Elona. ae 


which, because of Eq. [19.2] and Stokes’s theorem, becomes 


p{E+ 5x] ds = Bras. 


The theorem is thus proved. 


Chapter 4. Rapidly Varying Fields 


24. THE MAXWELL EQUATIONS 


The equations which we have applied to quasi-static fields 
contain an inconsistency. The equations for vacuum were 


4 
curl H= —i, divH=0, 
10H ; 
curl E=—~—, div E = 4no. 
From these it follows that 
divi=0. 


This, however, is not exactly true since, in general, we also 
have the continuity equation 


oO age —0. 
ot 


(Surface charge densities can always be considered as limit- 
ing cases of volume charge densities.) The situation can 
easily be corrected by adding an additional term. From 
the fourth equation and the continuity equation it follows 
that 
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Since div curl H =0, our equations become consistent if i is 
replaced by i+(1/472)(OE/ot). . Thus, 


The additional term (1/c)(OE/0t), which is called the dis- 
placement current, is due to Maxwell. This term vanishes 
for steady-state fields. We now have the following com- 
pletely self-consistent set of equations: 


cul E=—=H, [24.1] 
4 - 

curl H = = oe = E, [24.2] 

divH=0, [24.3] 

divE = 4%. - [24.4] 


These are the famous Maawell equations for vacuum. 

The continuity equation now appears as a consequence 
of the Maxwell equations, since from Eqs. [24.2] and [24.4] 
it follows that 


o+divi=0. [24.5] 


A somewhat analogous result holds for surface charge den- 


n, 
n, 
Figure 24.1 


sities. According to Eq. [7.3], 
4nw = E,,, ae ; 
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The conservation law for electrical charge requires that 
O+4, +4,,=0. 

From these two equations it follows that 
Aze(tn, + in,) + En, + Ep, = 0. 


That is, the normal component of 42i+E is continuous 
across a boundary. 


Molecular viewpoint. The above equations are correct for 
the field intensities in vacuum. They are, therefore, also 
valid for the microscopic fields e and h. Forming the spatial 
averages of e and h leads to the macroscopic fields E and B, 


e=E, h=B, 


for which the following equations hold: 


eurl E = — = B : [24.1'] 
curl B = =i ob -E, [24.2’] 
divB =0, [24.3] 
divE = 4x0. [24.4’] 


These are the Maxwell equations for material media. In these 
equations, @ denotes the total charge density (including 
polarization charges) and i is the total current density, 


i—i,tecurlM+P, [24.6] 


where i, is the conduction current density, ¢ curlM is the 
molecular current density, and P is the polarization current 
density. Since the polarization P is equal to the electric 
dipole moment per unit volume, 


re = (> ex) per unit volume ? 


P a OS ev) per unit volume 


then 
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is indeed a current density resulting from the variation of 
the dipole moment. 

The Maxwell equations [24.2’] and [24.4’] can ‘be written 
in terms of D and H, defined, respectively, by 


D=E-+ 42P [24.7] 
and 

H=B—4cM, [24.8] 
if we recall that @=0,—divP, as was shown in Section 7. 
This results in 


curl H = = ae =D ; (one 
div D = 479,. [24.4”] 


Just as Eq. [24.5] follows from Eqs. [24.2] and [24.4], the 
continuity equation for true (conduction) charges, 


@:+divi,=0, [24.57] 


follows from Eqs. [24.2"] and [24.4"]. 
To the above equations must be added the phenomeno- 
logical (semiempirical) equations: 


DSckE., [24.9] 
B= nH, [24.10] 
i, = cE. [24.11] 


The quantities ¢, uw, and o are not really constants; they 
are different for slowly and rapidly varying fields. For 
periodic fields 


F= F(x) ef ) 


where F represents any of the quantities D, E, B, H, or i, 
they are functions of the frequency: 


&(w), (a), a(w) . 


In the limit w—0, they converge to the ¢, w, and o for 
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static fields. The theoretical determination of the frequency 
dependence must be made with the aid of special conceptual 
models for the medium. In the final analysis, quantum 
theory must be invoked since the question of atomic struc- 
ture is involved. For arbitrary (nonperiodic) fields, e, y, 
and o are not even numbers but rather operators. The 
phenomenological equations [24.9], [24.10], and [24.11] 
are less general than the other equations. 

One consequence of the continuity equation [24.5’] and 
the phenomenological equations [24.9] and [24.11] is the 
following: 


divi, = odiv E =~ divD ser 
E 


: 420 
O.+ Wel 


Upon integrating, 


4no 


Qr=ore * 


Thus, a finite time is required for the true charge density 
in a conductor to disappear. This is called the relaxation 
time. 


25. ELECTROMAGNETIC WAVES IN VACUUM 


The existence of electromagnetic waves is a fundamental 
consequence of the Maxwell equations. In vacuum, 


e = 0 ’ CaS 0. 
From [24.1] it follows, by taking the divergence, that 


ka davi=0. 

ot 
Equation [24.3] is, therefore, not independent of Hq. [24.1], 
but it does, however, assert somewhat more. For periodic 


fields Eq. [24.3] is, indeed, even a consequence of Hq. [24.1]. 
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Analogously, Eq. [24.2], with i=0, leads to 


Cue 
a divE=0, 


and, for periodic fields, even to 
divE =o; [25.1] 


this latter result following in all generality from Eq. [24.4] 
when o=0. 

The additional term E/e in Eq. [24.2] has the physical 
consequence that electromagnetic waves can exist. That 
is, upon eliminating H from Eqs. [24.1] and 


curl H = ~E [25.2] 


by forming 


curl curl E = grad div E— V?E = — Ee ik 
Cc? 


one obtains, because of Eq. [25.1], the famous wave equation 
for E: 


flies 
V2E— a E=0. [25.3] 
Equations [24.1], [25.2], [24.3] and [25.1] remain correct 
upon the substitution 
E->H, H->—E. [25.4] 
Thus, it also is true that 


WH — H=0. [25.5] 


As an example of a solution of the wave equations [25.3] 
and [25.5], we consider plane waves: 


E = E,eth=-o H = Heh =- , [25.6] 
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(In these equations, the real parts should actually be con- 
sidered since only they have physical significance.) We 
have set 

O— 2m, [25.7] 


k=kn, [25.8] 


where y is the frequency and n is the wave normal; that 
is, the unit vector nm is perpendicular to the surfaces of 
constant phase, the so-called wave surfaces: (n-x) =con- 
stant. The phase also remains unchanged at a point a 
distance A=22/k away along the n direction. Thus, 


k= — [25.9] 


where A represents the wavelength. For solutions of the 
form given in [25.6], the differential operators are equivalent 
to certain multiplications: 


, 0 , 
Ag ~ik ’ at ~—tW, 
32 [25.10] 
Vi~—k, ace 
These substituted into Eq. [25.3] yield 
1 wo? 
(-# 45 w) H=0, k= ea [25.11] 


In order to see a constant phase, one must advance with 
the velocity w/k in the n direction, since then 


m= 4, k(n-x) = ot. 


Thus, 


a. 25.12 
aid [ ] 


is the phase velocity. The quantity c is contained in the 
differential equations and is independent of w. Furthermore, 
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it follows from Eqs. [24.3] and [25.1] that 
i(k:H)=0, itk-E)=0, 
so that 
n= nh 0. 


Thus both H and E are perpendicular to n; that is, elee- 
tromagnetic waves are transverse. 
Because of Eqs. [25.10], it follows from Eq. [25.2] that 


i(kx H) =—E, 


while from Eqs. [25.8] and [25.12] there results 
nxH=—E or E=Hxn. 


Thus E is perpendicular ton and to H. Hence, E-H=—0. 
Furthermore, 


|E|=|H| or E?=H?, [25 13] 


since n is a unit vector and H is perpendicular to n. Since 
the equations remain correct under the substitution given 
in Eq. [25.4], we have 


H =—(Exn)=nxE, 


and it is seen that E, H and n, in that order, form a right- 
handed orthogonal system. 


E 


Figure 25.1 
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There follows the relation 
Ex H = nF? = nH? = 3 (E?+ H?)n, [25.14] 


which we shall use later. 

These results obtained by Maxwell refer to electromagnetic 
waves. It is natural to consider light as being waves of this 
kind since 
1. the transversality, which is very difficult to explain in 
terms of a mechanistic aether concept, results quite natu- 
rally, and 
2. it turns out that c, which was defined as a pure electro- 
magnetic quantity and, as such, is determined by the mag- 
netic field due to a current, agrees with the empirically 
determined velocity of light. 

Light, therefore, simply signifies electromagnetic waves 
of certain special wavelengths (electromagnetic theory of 
light). Hertz and others later discovered electromagnetic 
waves of other wavelengths (radio waves, infrared). 


26. CONSERVATION OF ENERGY AND MOMENTUM 


We start with the Maxwell equations, [24.1], [24.2], [24.3], 
and [24.4]. 
According to Eqs. [17.5], the force per unit volume is 


il 
Thus, the work done per unit volume, per unit time, is 
Pe) [26.2] 


a. Conservation of energy 


Multiplying Eq. [24.2] by cE/4z yields 


see E-E—— E-cwlH= 0, 
4a 4n 
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while Eq. [24.1] multiplied by eH 47 gives 
1 . ¢ 
H-H+ -—H-eurlE=0. 
4x 4 


Adding these two equations and using the faet that 
divEx H = H-eurlE— E-eurl H 


results in the relation 
CW i 
A+ = 4+divS=0. [26.3] 
Cc 
Here, 
1 


Sa 


W= 


(B?+- H?) [26.4] 
is the energy density and 


Ss= 


c 

ExH 26.5 
4a [ } 
is the Poynting vector. Equation [26.3] is the law of eouser- 
vation of energy as it applies to field theery. Upen tntegra- 
tion, there results 


a ee — 
Jan +a ]was +O Sif = 0. [26.6] 


Here, the first term is the total work done per seeend by 
the system, and the second term the energy inerease of the 
system per second. Thus, the third term clearly must rep- 
resent the total energy flowing out of the bounding surtace 
of the system per second. Therefore, S has the siguitieance 
of an energy current density. 

If the bounding surfaces are absorptive or else are seo far 
away that the fields have net yet reached them, then the 
third integral vanishes. The work done by the system is 
then equal to its increase in energy. This is always the ease 
for static fields sinee then [S| at large distances falls off at 
least as fast as 1/r*. 


Sec. 26 | CONSERVATION OF ENERGY AND MOMENTUM 119 
For plane waves, according to Eq. [25.14], 
Sieve [26.7] 


This is very satisfying since it implies that the velocity with 
which the energy propagates is equal to the phase velocity 
of the waves. 


b. Conservation of momentum 


This consideration is similar to that of Section 18, al- 
though here there are additional terms. 

From Eq. [26.1] and the Maxwell equations [24.1] to [24.4], 
it follows that 


4nK = in| kav = ta| (E+ <ixH)av 


=|{E-aivB + curl Hx H— ~ Exh] dV. 
Because 
1. 
H-divH = 0 and (curl + 2H) xE=0, 
there results 


1 . 
0 =| {Heaiv + curlExE—~ExH} dv. 


With the abbreviations 


as =|{E-aivE + curlExE}dV, 


Jn=|{H-aivH + curlHx H}dV, 


we then have 


1d . 
4nK = Je +In— = 5; ExHdy. 


Now, it has been shown in Section 18 that 


Gq Uet In) =f Tod, 
aU 
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where : 
Te = Te + Te [26.8] 
is given in Eq. [18.7]. Thus, we obtain 
K= pT. df — a ae [26.9] 


The second term on the right is exactly zero for steady- 
state fields and is negligible for quasi-static fields. The 
differential form of this relation is 
OT, 1 08; 
k= >= = 
—“~ Ox,  ¢% ot 


(i=1,2,3). [26.10] 


The significance of Eqs. [26.9] and [26.10] is that momentum 
can be conserved only if a momentum as well as an energy 1s 
attributed to the electromagnetic field. If the total mechanical 
momentum is denoted by P,,, then 


== Kk [26.11] 
From this it follows that an electromagnetic momentum 
i 

a = /Sav . [26.12] 


must be ascribed to the field. This implies a momentum 
density of 


g= = S [26.13] 


per cubic centimeter. If the bounding surfaces of the region 
of integration in Eq. [26.9] are so located that at some definite 
time they have not yet been reached by the fields, then 


T.af=0, 


P [26.14] 
a (P+ Pom) =0, P,+ P,,.= constant, 


as shown by Poincaré and Hinstein. In general, however, 
P,, and P.,, are not separately constant. 
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c. Radiation pressure 


A plane wave is a mathematica] idealization: it is infinite 
in extent and has, therefore, infinite energy; its frequency 
is sharply defined. We can, however, allow the wavelength 
to vary within some arbitrarily small interval. This is suf- 
ficient to produce a finite wave train (theory of the Fourier 
integral). Such a wave train having very many humps but, 


GY 
(OLDGDOSS 


Figure 26.1 


nevertheless, being limited in extent is called a wave packet 
(Fig. 26.1). The energy and momentum are then given by 

i) E 
z=|way, Pa = 5 [SdV= 05, 
Ce? c 


according to Eq. [26.7]. Thus, a wave train whose dimensions 
are large as compared with its wavelength, propagating in a 
definite direction, has the momentum 


energy 


velocity of light 
This suffices for calculating the radiation pressure.! 


Simplest case. Total reflection at normal incidence 
(Fig. 26.2). The relation given by Eq. [26.9] is applied to 


1 The radiation pressure can also be calculated directly from Eq. {26.1}. 
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a surface at which the fields exist from a time t& until a 
time t,. The electromagnetic momentum changes by an 
amount 2#/e and a reaction is clearly exerted on the surface. 
The change in mechanical momentum is 


Apesee 
c 


If T denotes the length of time during which the wave train 
is in contact with the surface, then to sufficient accuracy, 


E=TFI, 


where F is the exposed area and J= S_ is the intensity of 
the radiation. Thus, 


AP er. 
yi 
where p=2 = [26.15] 


is the radiation pressure. 
It is easily seen that for a wave packet incident at an angle 3, 


iE 
p=—2 _ cos? @ . [26.16] 


Figure 26.3 
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For total absorption, 


[26.17] 


Figure 26.4 


These results cannot be obtained without the electro- 
magnetic theory of light. The essential fact that we have 
used is that the electromagnetic momentum per unit volume 
is S/c?. Radiation pressure was first experimentally detected 
by Lebedev. 

This association of momentum with the electromagnetic 
field is quite natural. Without it, conservation of energy 
and momentum at each instant of time cannot be obtained. 
If light is emitted and then absorbed, the conservation law 
for mechanical momentum is, indeed, once again valid after 
the absorption; however, it seems unnatural to eliminate 
the field since it is not apparent why more reality should be 
ascribed to the material particles than to the field. 

Thus, on the one hand, the aether (by aether we now mean 
all possible electromagnetic fields) is described independently 
without reference to mechanical quantities, and it is point- 
less to speak about a motion of the aether. On the other 
hand, the equations for a perfect vacuum (9=0, i=0) 
are only an idealization, since electromagnetic fields can be 
produced and detected only with the use of the mechanically 
describable particles that carry charge. A noteworthy 
duality is encountered here [A-1]. 
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27. ELECTROMAGNETIC WAVES IN MATERIAL MEDIA 


The differences between material media and the vacuum 
are attributed to magnetization and polarization. If we 
allow « and w to differ from unity, then the Maxwell equa- 
tions of [24.1'], [24.2’], [24.3’] and [24.4'] become, in the 
absence of true charges or currents (9,=i,=0), 


curlE =—"H, divii =, 
curlH=+-E, div — 0; 


The phase velocity of the waves is then 


c 


V eu 
For the index of refraction, defined by 
y 
Ve= ee {27 2| 
we obtain 
n= eu. [27.3] 


For plane waves, in analogy to Section 25, it is found that E, 
H, and n again form a right-handed orthogonal system. 
However, the ratio of the E and H amplitudes is no longer 
unity. Instead, 


i(kx E) = ik(nx E) = + "iol, 
[27.4] 


5 
\H|=|/= IE. 
be 
The last equation can be remembered most easily in the form 
eH? = wH’. [27.5] 
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With respect to Eq. [27.3], it may be remarked that para- 
and diamagnetism are always very weak, so that for non- 
ferromagnetic substances it is essentially true that »=1, 
n* = &(w), where lim n*(w) = é ra, (Static dielectric constant). 


This can be explained in terms of electron theory. For fer- 
romagnetic materials, uw is also essentially equal to 1 in 
the optical region of w since ferromagnetism does not occur 
at high frequencies. 


28. RADIATION OF ELECTROMAGNETIC WAVES 


In this section the generalization of the fields produced 
by static and quasi-static dipoles is described. 


a. The electromagnetic potentials. The inhomogeneous wave 
equation 


We start from the Maxwell equations in vacuum, 
[24.1], [24.2], [24.3], and [24.4]. From Sections 20 and 21 
we already know that Eqs. [24.1] and [24.3] are identically 
satisfied by the introduction of the potentials A and 9: 


H=curlaA, [28.1] 
ser 
E=— 34 = grad . [28.2] 
Introducing A and @ into Kq. [24.2] yields 
‘es 1 += 47, 
: 1.\_ o, 1 y_ 4a 
grad(divd += 9] VtA+ A aot 
because of Eq. [15.9]. Putting A and g into Eq. [24.4] 
gives 
— aaieed —V’%y = 420. 
C 
Here, it is useful to rewrite this equation slightly: 


10 thinks Noe 
mee Sia; ~¢\— V9 + = =4n0. 
C =, (div A i -*) es eal 
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It is now seen that it is convenient to require that 
div A + + ¢ = ( (the Lorentz condition), [28.3] 
; : 


rather than to demand that divA =0 as was done for the 
quasi-static case. Thus, one obtains the inhomogeneous wave 
equations for m and A: 


— 
NAD eg 0 = Ames [28.4] 


~vid+ 5 =i. [28.5] 


The three equations [28.3], [28.4], and [28.5] are consistent, 
since 
o+divi=0. [28.6] 


In vacuum (where go =0,i=0) these relations reduce to 
the homogeneous wave equations. 

The integration of Eqs. [28.3], [28.4], and [28.5] will be 
carried out in two different ways. 


b. Integration of the inhomogeneous wave equations. First 
method 


First of all, we assume that the fields are periodic in time. 
This results in an important simplification. Let the charge 
and current densities be 


O= Qe, = b= ie", [28.7] 
Then, 
(as Qe" ; A= A,e~ ; 


from which it follows that 


With 


—=k, [28.8] 
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the differential equations [28.3], [28.4], and [28.5] yield 


divA—ikp =0, [28.3'] 
Vp +k?o =—4n0, [28.4'] 
VA +A = — 2 28.5") 


These equations are easily integrated by the use of a 
method based on one employed in Section 6. There it is 
shown that the solution of 


— V2 = 4x0 
= "Q9°dVo 
es | Tro 
To obtain this result, Green’s theorem, 


- 7 oy op 
Jove y Vig) dt = (0 An vselar, 


is 


is used with y=1/r,, so that 
V7y=0 for PQ. 


For the present case one must set 


etre 
y= [28.9] 


L pg 


in Green’s theorem. This y satisfies the equation 
Veyt+ky=0 for r<0. [28.10] 


Proof: In polar coordinates Eq. [28.10] becomes 


Amare wiry] = 0, 


dr? 
which has the solution 


ry = eter, 
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(The reason that we have taken the +7 rather than the —7 
solution in Eq. [28.9] will be made clear later.) 

For the region of integration V, we again take the interior 
of a large sphere K about the point P. The point P itself 
is excluded by means of a small sphere k. Thus, because 
of [28.10] and [28.4’], 


e*" Pe 


gp Vip— pVip = Vip + key) — p(Vip + kp) = + 4709 
PQ 


This, substituted into Green’s theorem, gives 


tk? pg ikr tkr 
tn [oo P ae= (95, = - = =f af 
v 


Ype 
re eikr eikr oy 
+$(p 5, > ae 


Now, at the surface of the small sphere, 
0/on = —O/er, gar dQ, 
where dQ is the solid angle. Also, 
oe 1 


ee Nee Vet 


ke 


Thus, the integral over the small sphere becomes 
pao {0 ikrgg+r a eitr , 
or 
k 
Because 0/on=c/0R, the integral over the large sphere 


becomes 
-[ae le +R (e- ikg)| Erk | 


We now consider the limit in which the small sphere vanishes, 
so that r+0. Hence, 


[--s20,. 
k 
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As the large sphere becomes infinite, R-> oo and jim p= 0. 
In addition, we require that 


: op . )= 
lim R55 — ik ae [28.11] 


This is satisfied if one demands that in the asymptotic region, 
for large R, 


eikk 
oom 
Since 
PA id all describes outgoing waves 
aie describes incoming waves |’ 


the significance of Eq. [28.11] is that only outgoing waves 
are considered. (Had —ik been taken in the previous 
equations instead of +7k, then incoming waves would have 
been obtained.) Equation [28.11] is called the radiation con- 
dition (after Sommerfeld). This condition does not follow 
from the differential equations but is, instead, an important 
additional physical requirement. This condition must be 
imposed in order that the solution be unique. 

The integral over the large sphere now approaches zero 
as R-—->oo and there results 


ere 
Pp =[ee dV, C [28.12] 
Tq 
Using an analogous procedure for A, one obtains 
il e'*"Pe ; 
ae = fis dV. [28.13] 
c Tee 


The subsidiary Lorentz condition of Eq. [28.3'] is also 
satisfied because of the continuity equation [28.6] which, 
in our case, can be written as 


div i—iko= 0. 
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Employing the identity 
div (af) = fdiva + a-grad f 


twice, there results 


: il ; e'*Pe il ; a Fe 
div A = = [Veta grad, 7 =— efor ig: gradg a 
all . f era al : F ere 

=— = [8% divo (i — + ~ | aVelaive to) —_ 5 


The first integral vanishes when transformed into a surface 
integral by Gauss’s theorem and evaluated over the infinite 
sphere. Thus, 


ikr 


div A = i-k{ AV 00° = ike , 


Tre 
which is the Lorentz condition for periodic fields. 
Equation [28.3'] is, in addition, very useful for determin- 
ing » when A is known. If t—+t—r/e is substituted into 
Eq. [28.7], then, because of Eq. [28.8], we have the trans- 
formation = 
e7tmt 


aS e tettikr f 


The solutions of Eqs. [28.12] and [28.13] then can be ex- 


pressed as 
Pag GV 
gr(t) = i Go (1 “s) — [28.14] 
Ci) oy 
a all . lpq dV, 
Ant) = “fi (1- ‘| ee [28.15] 


This form of the solution is much more general than that of 
Eqs.[28.12] and[28.13]. It is not tied to a periodic (~e™”) 
time dependence. That is, if @ and i are linear superposi- 
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tions of periodic quantities (Fourier integrals), 


+c 
o(t) =| elope do , 


i(t) =|i(oye dw, 
then the corresponding solution is also a linear superposition: 
+0 
p(t) =|. emda , 


—o 


+a. 
A(t) =| Aloe dw. 


That is, Eqs. [28.14] and [28.15] are valid as long as 0 and i 
can be represented by a Fourier expansion in time. In the form 
given by Eqs. [28.14] and [28.15], my, and A, are called 
retarded potentials because of the appearance of t—r/c. 
They correspond to the radiation condition. The advanced 
potentials would be 


telt) = o0(t+ 2), 


Tre 


Arlt) = 2 fig(t + 722). 


Cc le 


These potentials represent incoming spherical waves. They 
are solutions of the same differential equations but are, 
however, not easily realizable in nature. In an infinite space, 
nature favors the first set of solutions. 


c. Discussion of the formulas for the periodic case 


A general consequence of the solutions of Eqs. [28.12] 
and [28.13] is the fact that a radiation zone is formed at great 
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distances from the current system (i.e., the region where @ 
and i are significantly different from zero). If the charges 
are restricted to a finite region of space, the static and quasi- 
static fields fall off as 1/r?. For the case of rapidly varying 
fields, not only the potentials but also the field intensities 
fall off as 1/r. 

For the following considerations it is important to dif- 
ferentiate between three different lengths: 
d = the linear extention of the current system; 
R = the distance from the field point Pto the current system; 
A = 2n/k= the wavelength. 


1. General discussion of the radiation zone. We now con- 
sider a field point P such that 


R>d and R>A 


(that is, kR>>1). These two conditions characterize the 


Figure 28.1 


radiation zone. Nothing at all is presupposed about the 
ratio d/A. 


If a point O within the current system is chosen as the 
origin of coordinates, then 
‘pa = R — Xo ’ 


where OP =R, X.=position vector of a source point Q, 
and n= x,/r is a unit vector in the P direction. In 
Eqs. [28.12] and [28.13] the denominator 7,, can be replaced 
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by & (this cannot be done in the exponent). Then, 
etkR 


aaa 


erkk 1 


R 


In the transition to the field intensities, 
oe 
a — curl, E=—-~A-—gradg, 


the differentiation results in two terms: in the differentiation 
of 1/R, the order of magnitude is decreased by a factor of 1/R; 
in the differentiation of the exponential function, multipli- 
cation by k results. Since kR>>1, the first term can be 
neglected and we have 


: «7 Xp : 
eyes —_ ~ eal, k * 
ay 1k Ox, ik R ikn 


Thus, 
H = curlA = tk(nx A) 
ikR 
ee 5g [nie eikin-20) dV . [28.16] 
Because 
div A = ik(n- A) 
we obtain 
g=n-A 


from Kq. [28.3’]. Thus, 

E = ik(A—n(n-A)) = tkA, 

tk eikR 
=a five ge Breede, [28.17] 

where 

i, =i—n(n:i). 
Since nxi=nxi,, it is seen that 
1. only the components of i perpendicular to n are effective in 

producing the fields E and H; 
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2. the same relations hold between E, H, and n as for the case 
of plane waves (Section 25): 


H=nxE, Hn=En=E-H=0, 


c Cc 
— 2 —_— — = 26 
|E|?}=|H, S=7ExH=—|E}-n 


Instead of Eq. [28.17] one can also write 


11 [(di, 
E=_=. | (Ft) a” [28.18] 


2. Dipole radiation in the radiation zone. In addition 
to the conditions R->>dand R> A, we now also require that 


A>d 


(that is, kd<1). This requirement is, for example, fre- 
quently satisfied for antennas. 

The exponential function e*"’* can be expanded in a 
power series. This corresponds to a decomposition of the 
radiation into multipoles (dipole radiation, quadrupole radia- 
tion, etc.). In the first approximation, one sets 


~ 


eT ik(n +x) =] . 


(That is, the retardation within the current system is ne- 
glected.) It then follows that 


ikR 
E = ik “fiay, 
‘ [28.19] 
erkR 1 
H= ik cnx fiay 
Re 


This is called dipole radiation. In this case, the dependence 
of the field intensities on n can be given explicitly. Let 


7 =fiav [28.20] 
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Then, 


ik eikR 
B= pds n= nxE. [28.21] 


For the case of a linear oscillation parallel to the z axis, 
lj,|=[j|sind. 
The Poynting vector (S = (¢/42)|E|?n) is then 


ce ke ue 
"dn oR? 


1 1/aj . 
= —<—<__—. —_ | — 2 
4nR? a(z)., we. 


sin? ? sin? (kKR— wt+ «) 


The energy radiated per second is obtained by integrating 
over the surface of a sphere: 


2 15fdj\* 
i dj = 3 3 (zy [28.22] 
since 
Al 2 1 1 
= = [eins == one 2 =a 
df hdao, 4, | 32 ba2 3) in | °° bdQ 3° 


The case of a radiating system consisting of a point charge 
is contained in the above if the charge oscillates as a unit: 


i= ov, j=fiav=ev. 
Thus, for an oscillating point charge, 


dz 2¢. 
[881-3 ga De [28.23] 


if A is large as compared with the linear dimensions of the 
region in which the charge moves. Naturally, A must also 
be large as compared with the dimensions of the charge: 
A>a. For an electron, a~e?/mc?~10-** cm. 
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It is often useful to introduce a polarization vector P: 


oP 


i=—, e=—divP. [28.24] 
ot 
For periodic processes, then, P=ii/w and 
etkR 1 oP, erkR 1 ail, 
ES -Fel ae meme: 
where 
n=|Par, j= [28.26] 


Nonlinear oscillations. Let II, and II, be two orthogonal 
linear components of the oscillation. Then, 


alia (1, = AW xn, 
erkk |] 
E, Felis then 


The Poynting vector is then not additive: 
c 
S=S,+S,+ qt et Exh}. 


However, the cross terms vanish when integrated over all 
directions. That is, their component in the n direction is 
c 


7 (n- E,x H,+ n:E,x H,) 


Cc 
a (E,-H,xn+ E,:H, xn) 
c 


2(E,:E,) = constant [II, — n(n-II,)]-[II, — n(n-Tl,)] 


Aa 
= constant [ (II, -II,)— (n-II,)(n-II,)], 
where constant implies independent of n. We now use the 


provision that II, is perpendicular to II,. With II, parallel 
to the a direction and II, parallel to the y direction, our 
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expression becomes 


constant x n,n, . 


This vanishes when integrated over the surface of a sphere. 
We thus have the important relation 


ps. af =o 81 af + 8x af, [28.27] 


if the oscillations are mutually perpendicular. 

Special case: circular oscillations. This case is important in 
optics. Here, linearly polarized light is emitted along the x 
and y directions, and circularly polarized light along the z 


y 


be 


/ 


Figure 23.2 


direction. In any other abitrary direction, the light is 
elliptically polarized. 


3. The Hertz vector. In the radiation zone, we had the 
condition kR>>1 (or R>A). We will no longer make this 
restriction. Instead, we now assume that 


a<i, ah, 


for arbitrary kR. This type of radiation is also called dipole 
radiation, although now we will no longer limit our consid- 
erations to the radiation zone. 

In Eas. [28.12] and [28.13] we can once again replace rp, 
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by an average value R. Then, 


1 etkk P 1 eth 


In order to show the connection with the usual notation, 
we set 


dll 
f= Ap = iwll ; 
so that 
—j=—ikil 
and 
Ls etkRk 
A=— tk — 1M he [28.28] 


Because of [28.3'], we obtain 


ikR 
o= —aiv (5 n) — — II-grad 


erkh 


R ? 


[28.29] 


since II is independent of position. 

The field intensities are now calculated without neglect- 
ing powers of 1/kR. It can already be seen that this solution 
is related to a dipole, since for k =0, m becomes the electro- 
static potential for a dipole of moment II. If, however, 
the dipole oscillates, additional terms also contribute signif- 
icantly. We have 


erkR 
H = curlA =—ik grad = x ; 
Because 


d 
grad j(k) = on, 


where n=x/R, then 


A d etkk . etkR s it 
H =~ ik () (ax = ik (ix ) mx. 
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Thus, 
erkk 


ae: 


(w+ 5 5 nxI1. [28.30] 


The first term has already been obtained in the radiation 
zone. The second term is large at small distances. Then » 
becomes 


d /eike : 
°-—agl@ (n- 1) = (- oe etkR xT]. 


Calculating grad yields 


tkR ike il y 
grady = 5, (-5 +7) + (i+ ots =) n(n) 


Thus, 


tkR 
E = ik-A—grady = — i T—gradg 


+(- ails zm" n(n) [28.31] 


This solution for H and E is known under the name of the 
Hertz oscillating dipole. 

The possibility of passing to the limit d>0 exists. Our 
equations are then solutions of the homogeneous equations 
(o=0, i=0) everywhere except at the origin, at which 
point a singularity exists. 

For the case where R<A(kR <1), one has 


li 
Ee R {3n(n-TI)—II}. 
This is the field of an electrostatic dipole of moment II. (This 


can be immediately seen by considering the z axis to lie 
along the n direction.) 
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For H there results 


tk al a 1ljxn 
a a + a 
This is simply the Biot-Savart field corresponding to the 


current j. 
In the radiation zone where R>>A(kR>>1), we obtain 
the previous results: 


etkk 
H=k R nxII, 
etkk 
E = — {H—n(n-T}. 


d. Integration of the inhomogeneous wave equations. Second 
method 


Almost all of optics can be treated with the method of 
periodic time dependence. However, we will demonstrate 
still another method for integrating the inhomogeneous 
wave equation 


V29 — — =~ = —Ano(x, t). [28.32] 


In this second method, no decomposition into Fourier 
integrals or Fourier series in time is made. 
A new variable t is introduced: 


i—t——, [28.33] 


where r=|x,—x,|. We also employ polar coordinates 
r, 9, y. Thus, 


r 
e(x;t) = o(r BY; 7-2) 


We denote 
OF (oF : ; ; 
a (=) (partial differentiation at constant t) , 
t 
28.34 
dF /oFr ee ‘ee 
apa (partial differentiation at constant 7) . 
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Then, 
OF adaF 10F ) 


6) 
ordre a’ & ar’ 
The divergence of a vector A becomes 


div,A=—A = (r?A,) + (derivatives with respect todand y) 


at fixed t, 


div, A= = 2 (r?A,)+ (derivatives with respect to @ and y) 


at fixed ¢. 
According to the above relations, 


: ; 1A, 
div, A = div, 4 —~ oan : 


Application to A= (1/r) V,y, where V, represents the gra- 
dient at fixed t, yields 


rt? 52 Or re dtor’ 


During integration, we shall require derivatives at fixed 1; 
thus we replace 0/or by d/dr. Then, 


wo (Uve)atvtp 10 118 18 dp_1 ay 
div, (=V.9) = ‘Pia dr reat reotdr re? ot 


deo: iL cl 10 
(Mem sae) mae (Pt Saf 2998 
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We will integrate this equation (at fixed t) over a volume 
from which the point P at r=0 is excluded by means of 


Figure 28.3 


a small sphere K (Fig. 28.3). Thus, 


op Oo(7, O, p; t — 1/C) 
[algae <=) oe 
FIK 
1 dp 1 op 
(e+e? 3) 22 +\(p ae nee 


where dQ is the element of solid angle and dV, = r?drdQ. 
In the limit K = the integral over the small sphere on 
the left-hand side vanishes, while that on the right-hand 
side gives 4ap,. With the aid of the geometrical relation 
dQ = (1/r?) cos (n-r)df, which is valid for any arbitrary 
surface F, there results the important formula 


we | gel”, 9, a t—#/¢) ay 


1 6) 10 
+2 ay |: (32) + (2 +2, cos (n OE [28.36] 


In this relation, t is fixed and r=7,,. 
This equation is also valid as a solution of the homogeneous 
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wave equation 


In this case, 


ee af {- (2) +(2+5 3) cos (ny r)} [28.37] 


If the surface F' is a sphere, then cos (n, r)=1, df= r?dQ, 
(0/On), = 0/dr, and there follows 


Any, =| dQ {¢ +r (2 tLe say: [28.38] 
FP 
For this case, 
lim gy,=0. 


To 


If, in addition, the radiation condition 


, op , 1 op 
limr & mn +B) 0 [28.39] 


is satisfied, then the surface F' can be taken at infinity so 
that the surface integral vanishes. 

For p=Me ‘, k=w/e, the previous relations are once 
again obtained. 


29. THE FIELD OF A POINT CHARGE IN UNIFORM MOTION 


The consideration of this case provides preparation for 
the special theory of relativity. 
We start with the equations 


diva += § 07 
ieee 
Vip — = =—4n0, 


VA 5 A=. 
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Let the charge 9 move with velocity v in the « direction so 
that #,=vt. We seek a solution of the Maxwell equations 
which depends only upon the position of the field point 
relative to the position of the charge, (« — vt, y, z). The 
time dependence, therefore, should appear only in the form 
z—vt. Then, 


From 


it follows that 


A,=A,=0, —=~*—-=0 
This is satisfed by the expression 
= (v/e)p . 
Only the equation for 9g, 


(ie ee Op dp _ 


ant gy: ee Ano , [29.1] 


must still be satisfied. The factor (1—v?/c?) implies a 
contraction. In order to obtain the same result as for the 
static case, we introduce 


,_ &—vwt 
The equation for y then becomes 
Op ep oe Op 
aa's T Sys ba ea —4n0, [29.3] 


which has the solution 


g=-—. [29.4] 
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Here, however, 
(3 css (x ae vt)? 


= Bet + y? + 22, [29.5] 


For the electric field it follows that 


1: 
E =— 7A grado ; 
vdp op v*\ op 
gee ES ay fs ea tek 
Cor ce (1 5] Oa’ 
Op oy 
E, Oy’ a 02’ 
et Vat 
7 
=" (ee [29.6] 
Co r 
@ 2 
errs 


Therefore E is in the direction of r=(#—vt, y, z) and not 
of r’=(a',y,2). The equipotential surfaces 7’= constant 
are flattened ellipsoids of rotation (Heaviside ellipsoids). 
The vector E is not normal to these equipotential surfaces. 

For the magnetic field it follows that 


With v= (v, 0,0), this can be written as 
H— “ voxE [29.7] 


and 
mie [29.8] 
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These relations are the precursors of a much more general 
result: The Maxwell equations are invariant under a group 
of transformations (the Lorentz group) which transform 
between coordinate systems in uniform motion and have 
the property that 


wl? y-+-2t— ¢ ogre — e+ y? + 2F— ce? + 2 2 


For a special choice of coordinates, such a transforma- 
tion has the form 


le %— vt lay g=—z j= t— (v/c*)a 
V1I—v2]o2’ V1—v/0 


The field intensities also transform in such a manner that 
the Maxwell equations remain unchanged. The laws of 
nature are thus the same in all reference systems. What 
is new is the fact that the time also undergoes a transforma- 
tion. In Newtonian and Galilean physics one has 


v—~—o, _=Sy, ear, tet 
Here, on the contrary, two events at different space points 
which are simultaneous in one reference system are not 


simultaneous in another. This is very closely related to 
the finite propagation velocity of electromagnetic fields. 


30. RADIATION DAMPING 
We have seen (Eq. [28.23]) that the energy radiated per 
second by an oscillating point charge is given by 


nee oe [30.1] 


According to the law of conservation of energy, the field 
must exert a reactive force on the charge. This force can 
be guessed at by considering the fact that 
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Since, however, the energy fluctuates, one need only demand 
that this be true for the time-averaged quantities: 


The first term on the right vanishes, since it is a total 
derivative. Thus, the radiation reaction (the radiation 
damping force) is obtained: 


eee ee [30.2] 


This is only a heuristic consideration. However, this force 
can be derived directly from the reaction of the field on an 
electron. For this purpose, the charge must be assumed to 
have a finite extent. It is found that Eq. [30.2] is correct if 
terms of the order of (v/c)? are neglected. 


Calculation of the self-force on a charge of finite extent [A-4] 


The electron will be represented by a charge distribution 
in motion, (9, 4), regarding which certain special assumptions 
will be made. 

The field produced at an arbitrary point P by the electron 
is given by 


te 
E, = —7, Ap— gradege [30.3] 
i, — curl,A,, [30.4] 
where 

Tpo\ AV, 
galt) =f oo(t—"), [30.5] 

PQ 

all 6 Ypa dV, 

= ps | a 30.6 
Ant) == fig (t— 22) [30.6] 


We are interested in the field existing in the region of the 
particle itself. Since this region is small, the retardation 
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will not play a major role. We expand in powers of —r/e: 


ti n=0 ¢ 
oud 1 es 
ae ae ot” | Tre 2a Ve, 
121 "8" fy 


In the expression for g, the n=0 term is simply the 
Coulomb potential. The n=1 term vanishes since i 0, dV, 
=e=constant and d(constant)/ot=0. In the remainder 
of the series we replace n by n+2. Thus, 


1\2t2 Ont2 es 
raf ean 3 otan(—) eh aren, 


Here, 0/ot operates only on g,. Using the continuity equa- 
tion 
2G¢ 
ot 
the identity of Eq. [4.10], Gauss’s law, and the fact that 
under the integral sign, 


= — divgig, 


a Io -gradg i ig ‘grad, 9 
we get 


| a a o 2+1 
Pr = Vo— a2, nee == ar ig? gradp rng adV,. 


The electric field is then found from Hq. [28.2]: 


0g NS Ota 
- a,| £ av,—= ¥ a) eal 
gra | Pap Sl ) ott! 


ae n 
“fle ign, aaa ore p(igogtads ne) adv... [380.7] 


In order to calculate the force exerted by the field on the 
particle we make the following assumptions: 
1. The charge distribution is rigid and moves only with 
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translational motion (although not uniformly). That is, 


i(x, t) = o(x, t) v(t), [30.8] 
where v(t) is independent of x. 

2. The reference system is chosen so that v/e<1. We keep 
only terms of zero order in v/e. This places no restrictive 
assumptions on the magnitude of the acceleration or on 
higher time derivatives of v. Indeed, our derivation is exact 
in the instantaneous rest system of the electron. Further- 
more, terms of first order in v/e vanish so that the first 
neglected terms are of order (v/c)?. 

According to Eqs. [26.1] and [30.8], the force per unit 
volume is 


1 
k= o(E+exH| . 
As the total force on the electron there remains 
K =) rE, dV, ’ 


since the magnetic Lorentz force is of order (v/e)?. If the 
expansion of Eq. [30.7] is substituted into the above expres- 
sion for K, it is seen that the Coulomb term drops out since 


Xp— & 
| foro 73 dV,dV, = 0 
PQ 


(as can be seen by interchanging the variables of inte- 
gration). There remains 


K=3.K, [30.9] 


n=0 


where 


“1\nt2 ont 
K™ = — (- ;) [ovate yeti 


1 1 , 
x [ee [one (w+ 2)! grad, (o-grady rf) aV,. [30.10] 


The principal terms are K® and K‘”. 
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In the calculation of K®, a term with 0@,/ot appears. 
Because of Eqs. [24.5] and [30.8], this is equal to —div,(o,v). 
This term, 


+ — a = | [ee divg (Ov) jv {v= — —5arad, (v-grad, r10) aV,dV,, 


can be transformed into 


— 5 | [ercer: grad, {v= — 5 erad, (wv: gradp,)| dV,dV, 


by using the integral theorem of Eq. [8.7]. Since the dif- 
ferential operators affect only r,,, grad, can be replaced 
by — grad, and vice versa, and there results 


; je 
e=- feece»-erad, {v= — 5 ~pradg (v: gradgrre) ay Su, . 
PQ 


The same result—although with the opposite sign—would 
be obtained by interchanging the variables of integration. 
This shows that the term must be zero. There remains 
for KK, 


Ko — — | fee Oo {o = ae 1 scl (o “grader av avy. 
Pe 


This integral can~-be simplified by assuming a spherically 
symmetric charge distribution. The integral of 
O*F'p9 it = (Xp, — Tq,) (Lp, — Lo,) 


= Ver 


a = 2 
OLp, Op, Tp Tre 


vanishes for iAk, and 0*r,,/ox, can be replaced by 
V?r/3 = 2/3r. Hence, 
: oe. 
grad p (v-gradprpg) > =v — 
3 1p 
and 


K® = — yo, [30.11] 
where 


21 
bas af [ee dV, dVo. [30.12] 
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Since 


Alen dV,dV, = Eu 
2 Pe 


represents the electrostatic energy (see Eq. [4.5]), can 
also be written as 


H=s—. [30.13] 


Because K has the character of an inertial force, u can be 
interpreted as an electromagnetic mass. If the charge is con- 
sidered to be localized within a sphere of radius a, then E,, 
is of the order of magnitude of e?/a, where e is the electron 
charge (cf. Section 5). 

For K® there results 


i oe 1 
Ko: a Gaels Pat {» e grad, (v -grad,r*,) dVo. 


Since 
grad, (v: gradprz,) = 2v 


and 
Joear, =| eua¥e= é, 
therefore 
homeo [30.14] 
3 c3 


is the radiation damping force. It is completely independent 
of the charge distribution. (In particular, a—0 is permis- 
sible.) 

The higher terms in the expansion of Eq. [30.9] are pro- 
portional to powers of a: 


K®n~ a - K® ~ a. 


These can be neglected if a is sufficiently small. 
Thus, in addition to the damping force given by Eq. [30.14], 
still another term, that of Eq. [30.11], is obtained. However, 
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if the equation of motion of the electron, 


,=—K 4 Wee eesre ee 
my = A =— pw ora 
is written, then it is seen that m and w, the mechanical and 
electromagnetic masses, respectively, cannot be separated. 
since this equation can also be written as 
eng eager 30.15) 
ye = ae as [30.15 | 
Hence, m-+yu is simply the total mass of the electron. 
For a point electron (a=0), the electromagnetic mass yp 
is infinite. 
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Appendix. Comments by the Editor 


[A-1] (pp. 2, 3, 6, 123). These passages, although not 
written by Pauli himself, reflect Pauli’s deepest convictions 
about the fundamental problems of theoretical physics, 
convictions that he repeatedly expressed during all his 
scientific life. This observation is substantiated by the 
following references to Pauli’s works (volume and page 
numbers in square brackets refer to Collected Scientific 
Papers by Wolfgang Pauli, edited by R. Kronig and V. F. 
Weisskopf, John Wiley & Sons, Inc. (New York, 1964)): 
Verhandl. Deutsche Physik. Ges. 21, 742-750 (1919) [vol. 2, 

pp. 8-9]. 

Encyclopddie der mathematischen Wiss., vol. 5, Part. 2, 

B. G. Teubner (Leipzig, 1921), pp. 539-775 [vol. 1, p. 237]. 
Scientia 59, 65-76 (1936) [vol. 2, pp. 747-748]. 

Prix Nobel 1946, pp. 131-147 (Stockholm, 1948) [vol. 2, 

p. 1095]. 

Helv. Phys. Acta, Suppl. 4, 261-267 (1956) [vol. 2, pp. 1304- 

1306]. 

Universitas 13, 593-598 (1958) [vol. 2, p. 1367]. 


[A-2] (pp. 29, 37, 41, 55). The general electro- and mag- 
netostatic property of matter is to be dielectric and diamag- 
netic, e>1 and uw<1. This is reflected by temperature 
independent « and mw (except for the diamagnetism of 
semiconductors), that is, e and yw are not of statistical origin. 
This case corresponds to the second analogy of p. 57; and 
e—1>0 and (1/u)—1> 0 express the fact that work has 
to be supplied to polarize the medium. 
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In the case of molecular or atomic (but nonpermanent) 
electric or magnetic moments, the first analogy of p. 57 
applies; e—1>0 and ~—1>0 are due to alignment in 
the applied field and again express the fact that work has 
to be supplied to polarize the medium. Here e and yw are 
of statistical origin and therefore temperature dependent. 

In the case of a permanent electric or magnetic moment, 
P or M is a nonlinear function of E (ferroelectricity) or 
H (ferro-, ferri-, antiferromagnetism), respectively. 


[A-3] (pp. 64, 67). The boundary-value problem for a 
steady-state conductor, og/on =0 on the lateral surface, 
yY = 9, = constant and »y = gy, = constant on the basal sur- 
faces, has a unique solution g. The proof is exactly the 
same as the one for the electrostatic problem (p. 41). 

For cylindrical symmetry, y is independent of azimuthal 
angle (also called g on pp. 65, 67), and since o = constant 
one verifies by inspection that g = — E-z, E = constant, 
is a solution. Hence i= constant. 


[A-4] (pp. 147-152). In an interesting recent contribution 
to the theory, Rohrlich remarked that the conventional 
calculation described here and elsewhere in the literature is 
not covariant and therefore [30.13] is not the correct ex- 
pression for the electromagnetic mass (F. Rohrlich, Classical 
Charged Particles, Addison-Wesley Publishing Co. (Reading, 
Mass., 1965)). The correct procedure yields, according to 
Robhrlich (his Eq. (6-21)), 


Prat = MO , 
where P%,, is the part of the mechanical momentum 


4-vector of the particle due to self-interaction, v* the 
velocity 4-vector, and 


_ Ba 
alr 


This is much more satisfactory than [30.13] since it satisfies 
the Einstein mass-energy relation. 
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Ampere, “fundamental law” of, 81 

Amperean currents, 74 

Analytic function, 51 

Asymptotic region, 129 


Biot-Savart field, 140 
Biot-Savart law, 7], 86 
Boltzmann, L., 1 
Boundary conditions, 

at conductors, 17, 18, 40 

at insulators, 40 

Boundary value problem, 156 


Capacitance, 96 
Cauchy-Riemann equation, 51 
Charge, 

conduction (true), 30, 112-113 
elementary, 2 

point, 4, 6-11, 18, 135, 143 
polarization, 30, 32, 33, 36, 112 
surface, 16-18, 19, 110 

test, 2, 3 

total, 111 

volume, 14-15, 19 
Charge density, 

surface, 16, 30-32 53 

volume, 15, 30, 32-33 
Circuits, 92—103 


Clausius-Mosotti relation, 39, 50 
Condenser (capacitor), 28, 96, 97, 
103 

Conductivity, 60, 61 
Conductor (metal), 17-18, 40-42, 

45, 47, 48, 49, 56, 59-81, 93, 95, 

103, 106-108 

Conservation, 

of angular momentum, 85 

of charge, 5 

of energy, 7, 117-119 

of momentum, 119, 120, 123 
Conservative field, 7 
Constitutive equation, 55, 112 
Continuity equation, 60, 109, 129 
Coordinates, 

curvilinear (orthogonal), 43-44 
cylindrical, 37-38 

polar, 24-26 

Coulomb (unit), 58, 127 
Coulomb’s law, 4—5, 28 

Current, 

conduction, 75, 76, 111-113 
convection, 75, 76 

definition of, 59 

displacement, 110 

molecular, 74—75, 11] 
polarization, 111 

steady-state, 59-81 

total, 11] 

Current density, 

of charged particles, 61 
definition of, 59 

surface, 74-75 


158 


Current density (continued) 
volume, 74-75 
Current distribution, 64 


Diamagnetism, 55, 125, 155 
Dielectric, 28-40, 46, 50, 155 
Dielectric constant, 29, 37, 39-40, 
Son lZ5 
Dipole, 
electrostatic, field of, 34, 139 
Hertz oscillating, 139. See also 
Hertz vector 
Dipole moment, electric, 33, 36, 39 
Ill 
Dipole radiation, 124, 137 
Discontinuity, 
of E at charged surfaces, 17 
of induction, 76 
of parallel component of D, 33 
Displacement current, 110 
Displacement vector, 32, 36 


Effective field, 36 
Einstein, A., 120 
Einstein mass-energy relation, 156 
Electric and magnetic quantities, 
analogy of, 57, 77, 155 
Electric dipole, 33 
Electric dipole moment, 33, 36, 39, 
111 e 

Electric field, 

definition of, 5 

of dipole, 34 

macroscopic, 37, 111 
microscopic, 37, 111 
Electrodes, 64-65, 94 
Electromagnetic mass, 151, 152, 156 
Electromagnetic waves, 

in media, 124-125 

in vacuum, 113, 117 
Electromotive force, 95, 99—100 
Electron theory, 2, 33 
Electrostatic potential, 8, 40—42 
Energy, 

of current system, 90—92 

of electric field, 18, 58 

radiated, 135 
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Energy current density, 118 
Energy density, 
of electric field, 21, 43, 86 
of electromagnetic field, 118 
of equipotential surfaces, 18, 145 
of magnetic field, 86 


Faraday, M., 1 
Faraday’s law of induction, 88-90 
Ferromagnetism, 55, 125, 156 
Field concept, 1, 7 
Force, 
acting on a charge, 80 
between two point charges, 4 
on a conductor, 78 
electric, 79, 80 
electromotive, 95, 99-100 
lines of, 1 
Lorentz, 79 
magnetic, 79-81 
ponderomotive, 77 
Forces, superposition of, 5 
Frequency, 115 
Function theory, 51 
“Fundamental laws,” 80-81 


Gauss’s law, 9 

Gauss’s theorem, 10, 11—14, 35, 82 

Gilbert, W., 54 

Giorgi (units), 58 

Grassman, “fundamental law” of, 
80 

Green’s theorem, 27, 127 


Heat conduction, 104 
Heaviside (units), 58, 86 
Heaviside ellipsoid, 145 
Henry (unit), 94 

Hertz, H., 1, 117 

Hertz vector, 137, 139 


Images, method of, 45-51 
Index of refraction, 124 
Inductance, 

mutual, 93 

self-, 92 

self-, for solenoid, 93 
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Induction, law of, 88-90 
for moving conductors, 106-108 
Insulator, 29, 34, 40-42, 48 


Joule heat, 63, 65 
Joule’s law, 63-65 


Laplace equation, 11, 51 
Laplacian in curvilinear 
coordinates, 43 

Larmor, J., 2 

Lebedev, P. N., 123 
Light, 

electromagnetic theory of, 117, 123 

polarized, 137 

velocity of, 117 
Liouville’s theorem, 51 
Lorentz, H. A., 2, 79 
Lorentz condition, 126 
Lorentz forte, 79, 149 
Lorentz group, 146 


Magnet, 106, 107 
Magnetic and electric quantities, 
analogy of, 57, 77, 155 
Magnetic charge, 54 
Magnetic dipole, 54 
Magnetic field, 
definition of, 54 
of dipole, 56 
macroscopic, 111 
microscopic, 111 
Magnetic induction, 54, 76 
Magnetic moment, 73 
Magnetic monopoles, 54 
Magnetic permeability, 55, 58 
Magnetic pole, 77 
Magnetic ponderomotive forces, 78 
Magnetic surface density, 55 
Magnetic susceptibility, 87 
Magnetic volume density, 55 
Magnetization, 54, 74 
Maxwell, J. C., 1-3 
Maxwell equations, 
in absence of charges or 
currents, 124 
for material media, 111 
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for vacuum, 110 

Maxwell stress tensor, 85 

Metal. See Conductor (metal) 

Method of images. See Images, 
method of 

Momentum, electromagnetic, 120, 
123 

Multipoles, 134 


Oersted’s law, 65, 86 
Ohm (unit), 58, 62 
Ohm’s law, 60, 106 
Oscillations, 
aperiodic (overdamped), 112 
damped, 101, 102 
forced, 100, 103 
free, 101 


Paramagnetism, 55, 125 
Periodic time dependence, 99, 104, 
112, 126-131 

Permeability, magnetic, 55, 87 
Permutation symbol, 12—14 
Phase velocity, 115, 119, 124 
Plane waves, 114, 119, 124, 134 
Poincaré, H., 120 

Point charge, field of, 6, 9, 136, 143 
Poisson equation, 15, 24 
Polarizability, molecular, 36, 39 
Polarization, definition of, 25 
Polarization charge, 30 
Polarization current, 111 
Polarization vector, 36, 136 
Polarized light, 

circularly, 137 

elliptically, 137 

linearly, 137 

Ponderomotive forces, 77—78 
Potential, 

advanced, 131 

of dipole, 34 

electrostatic, 8, 40—42 

retarded, 131 

Potential energy, 

for point charges, 18 

for volume charges, 19 
Potential flow, 64 
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Poynting vector. See Energy ampere, 58, 62, 86 
current density coulomb, 58 
electromagnetic (emu), 59, 62, 65 
Quasi-static fields, 88, 132 Giorgi, 58 
Heaviside, 58, 86 
Radiation, dipole, 134, 137 henry, 94 
Radiation condition, 129, 143 ohm, 58, 62 
Radiation damping force, 147, 151 rationalized mks, 58, 87 
Radiation pressure, 12]—123 volt, 58, 62 
Radiation reaction, 147 watt, 63 
Radiation zone, 132, 137 weber, 87 
Relaxation time, 113 
Relativity, theory of, 2, 143 Variation principle, 45 
Resistance, Variational problem, 42, 43 
definition of, 62 Vector potential, 
specific, 61 definition of, 55 
Resonance, 100 of field of dipole, 56 
Retarded potential, 131 Velocity of light, 117 
Rohrlich, F., 156 Volt (unit), 58, 62 
Rotations, 
improper, 12 Watt (unit), 63 
proper, 12 Wave equation, 
homogeneous, 114 
Scalar, 14, 86 inhomogeneous, 126 
Self-force, 147 Wavelength, 115 
Self-inductance, 92, 93 Wave normal, 115 
Skin effect, 103-105 Wave packet, 121 
Solenoid, 67, 93 Waves, 
Sommerfeld, A., 129 _ incoming, 129 
Stokes’s theorem, 8, 11-14 outgoing, 129 
Stream lines, 64 plane, 114, 119, 124, 134 
Stress tensors, Maxwell, 85 spherical, 131 
Superposition of forces, 5 Weber (unit), 87 
Susceptibility, magnetic, 87 Work, 7, 89, 90, 155 
Tensor, 


antisymmetric (skew-symmetric), 
12-13, 66 
concept of, 81, 82 
Maxwell stress, 85 
symmetric, 82 
Torque, 85 
Transversality, 117 


Units, 
absolute electrostatic charge 
(esu), 57, 58, 60 
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